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ABSTRACT 


This dissertation looks into the proble® of conically 
cambered delta wings with subsonic leading edges in a steady 
supersonic stream. The superiority of subsonic over supersonic 
leading edges has long been recognised because of its favourafele 
lending edge suction which alleviates drag. This provided the 
motivation for obtaining analytical solutions for the indirect 
lifting problem by choosing a tAvlst distribution function 
which is a polynomial of even ordered terms of the conical 
variable. Thus, wing shapes cambered in a plane perpendicular 
to the free-st ream are easy to visualise and definite trends 
in the results may be established. 

Considerable attention is paid to the problem of obtain 
ing minimum drag shapes subject to a given lift and other 
equality and inequality constraints. In order to' avoid flow 
separation at the leading edges, the attachment lines were 
fixed there. 

A noted contribution of this dissertation is the use of 
mathematical programming techniques to formulate the minimum 
drag problem subjected to realistic constraints which may quife 
often be stated as inequality constraints. For programming 
methods, it is shown that orthogonal loadings , as defined by 
Graham, play an Important role in obtaining a simplified fortnu. 

lation. It has also been possible to generalise some drag 

' .■ * 

reduction theorems already known in the literature. 



Vi 

The olf— design conditior?, when Adir.e voytioes 

are observed if the wing is slender and the leading edges 
are subsonic, is studied theoretically, a hybrid method is 
proposed which uses Polhamus’ leading edge suction analogy 
and a modified Kuchemann-Squlre forraul ation and the problem 
Is solved as a minimum drag problem. Reasonably good agree- 

nent at low aspect ratios with available experimental data 
is observed. 
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CHAPTEB 1 


INTm>I5UCTI0N 

i,i General Retnarkg 

Delta wings bawe frequently found applications in 
tbe design of supersonic airplanes because of its ease in 
manufacture, high wing thickness near the root ohord even 
for small thickness-chord ratios and greater structural 
rigidity. Fortunately, it is also a planfom suitable for 
theoretical analysis aiu! has been widely studied « 

This dissertation studies the analysis of a delta 
wing, with subsonic leading edges, flying at supersonic Mach 
numbers and hawing a conical twist distribution, fhe super- 
iority of subsonic over supersonic leading edges for such 
wings has long been recognised because of Its favourable 

ft 

leading edge suction which reduces drag. Further, any delta 
wing must pass through a Mach number rai^e in supersonic 
flight for which the leading edge is subsonic. The problem, 
therefore, is of practical interest. 

Conical comber was chosen because, within this ffreBse- 
work, it is possible to obtain closed form solutions without 
requiring the wing to be slender. Both the direct and the 
Indirect prohleras have been attempted In the literature, w© 
have chosen to solve the indirict lifting problem by specify- 
ing the twist distribution In even powers of the conical 



coordinate 


Since wing shapes of mintnmiti drag are of considerable 
practical imf»ortance, studies have been made using corabina- 
tions of these twist distributions. In the literature miniraiun 
drag wings have generally been obtained for a given total 
lift and, sometiRses, for a given pitching juoment. These con- 
straints have always appeared in equality form. The procedure 
is extended to include both equality and inequality constraints 

by using the powerful techniques of mathematical progra^^'mlng . 

32 

Fn the process the fisefulness of orthogonal lo<»dings is shown. 

As a corollary to this work it was possible to generalise 
some of the existing, drag miriiraization theorems by relating 
the forv.'ard and reverse flows. Those theorems, it is felt, 
shoTild he of help in not only identifying niiniiauti! drag sluijtes 
but also in obtaining wings with low lift dependent drag. 

At oil-design conditions, if tlic wing is slender, 

witi! sharp leading edges which are subsonic, leading edge 

vortices give rise to additional lift which is a non-linear 

function ol tfie incidence. A now tiieorctical model which uses 

Polhanus* suction analogy'^ ^ and a modified version of igthe 
54 55 

K«chemann‘ , Squire' ‘ formulation is used and the problem sol- 
ved as a drag minimization problem to iiredtct the pressure 
distribution. AgreeiWeut with available experinental data is 
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1 . 2 General Review of Supersonic Vi'in/r. Theory 

Conventional airfoil shapes fortunately allow us to 
linearize the equations of motion and to simplify the boundary 
conditions. Insplte of this, the problem of a finite wing in 
STibsonic flow has been quite difficult. The case of supersonic 
flow has been qiiite encouraging and in several instances have 
given closed form solutions. 

As one moves from the subsonic to the supersonic case, 
a distinct change in the character of the disturbance field 
takes place. The fore and aft symmetry of the subsonic flow 
disappears and the disturbances are limited to the region behind 
the wing bounded by Mach cones. Because of this asyr;metry in 
the flow all outgoing disturbances result in a drag, called the 
wave drag, so that even the tluckness of a wing or body creates 

t 

a drag in addition to the drag associated with lift. Thus the 
drag arising from the lift may be identified piirtiy as vortex 
drag and partly as wave drag, though both effects are due to 
the rearward tilt of the lift vector. The vortex drag, which 
in subsonic flow is also known as the induced drag, be cal- 

culated from considerations of the vorticity in the waive, and 
depends solely on the spanwise distribution of lift. i’he wave 
drag annoars due to the energy required to extend the wove 
systeni emanating from the wing to infinity. 

■^0 

Hayes*' shoived that at great distances behind the wing, 
the wave disturbance is concentrated in the vicinity of the Mach 




VORTEX WAKE 



MACH CONE 


FIG. I.l. DiSTURBANCE PRODUCED BY WING AT SUPERSONIC SPEEDS. 



5 


cone and docs not modify the purely subsonic field of dis- 
turbance in the vicinity of the vortex wake, which occup.ics 
a relatively small region near the centre of the cone. 

Figure 1.1. Ifonce tlic wave drag and the vortex drag may 1)e 
separated if coJculatlons are made on the basis of the nronien- 
tum of the flow at a large distance hehin<l tbe wing. Alter- 
nately, the sut'i of the wave and the vortex drag or the drag 
due to lift, may be calculated by integratinfr the product 
of the local pressure with the local angle of attack over the 
entire wing surface. 

Using the latter approach it can readily be shov/n 
that tlic drag thie to tfil cicncss is directly proportional to 
the square of the thickness-chord riitio , and the drag due to 
lift is i>ropor tionul to the square of the lift coefficient. 

In the linearized theory, it can he further shown that the 
drag due to thickness ami llft’are additive and hence the 
tw'o problems can be treated independently. 

Important distinctions arise between cases where the 
velocity component normal to an edge is subsonic or supersonic 
The edges are then known us subsonic or supersonic edges res- 
pectively. For a lifting wing with a subsonic Ic^ading edge, 
the flows on the upper and lowci* surfaces near the edrre are 
no longer inilepeiulent . The fluid flows around the edge with 
subsonic velocity, and one obtains a leading edge suction. 
l'urthernor<‘ , in the sector enclosed hetw'een the leading (uige 



and the limiting Mach line, the pressure restilting iron tin? 
upper and lower half spaces must be balanced. In the case 
of a subsonic trailing edge, the Kutta condition must be 
satisfied. 

In the case of a supersonic leading edge, the flow 
will be locally similar to an Ackeret flow, provided the 
edge is Shari'. There will be no flow around the edge and 
the surface of the airfoil divides the flo^v into itulepcnrient 
u"pcr and lower r€?gions. Tt is also no longer necessary to 
distinguish betwcon pressures arising from lift and those 
arising from thickness, 'fho pressure on each stirface being 
determined by the shape of tiiat surface alone. 

Tn the case of a supersonic trailing edge, since 
disturbances downstream of tlie trailing edge cannot propagate 
upstrean, the independent pressures on tlie u- (per ‘and lower 
surfaces will bo maintained upto a point so near the trailing 
edge that eq'ializatiou of the pi'cssnrc takes place throu;r.]i 
the bfumdary lay fir anmev/hat gradually. ilowevcr, in inviscid, 
linearized theory, it is c:ojisidered sufficient to let the 
lift fall to zer^> behind the trailing edge through an attached 
oblique shock. 

The case when both the leading and trailing edges 
are supersonic i s th c simplest to solve, since the flows on 
the top end bottom surfaces are com'vlctely independent and 
the lifting and thickness nroljlcns may be tr<'at(?d in the sane 
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fashion. 'I’his Is done hy rcMresentin." the surface hy a dls- 

trih»ition of stipersnnlc sources, u concc]»t first introduced 

by von Karnaii and yoore for axisy= metric 1>odies and later 

‘-*2 

extended to ivin^s by l’ucl';ctt‘“‘7 Tn the linearized theory it 
was shown that the local source strtfiif^th Is j»r<ti)Ortional to 
the local an,"le of attack or tlie str<?anmise sione. The cal- 
culation for the pressure distribution then reduces to integ- 
rating a double integral ratlier tlian solving an integral 
cfpuition. 

The 'probl em of partially subsonic leading edges was 

24 

elegantly solved, apparently independently, by Kvvard and 
25 

J'rasi lchlkova‘“ . The method requires at least some part 
of the leading edge to bn supersonic. Tn their method, to 
obtain the pressure at a point, it was necessary to integrate 
over the shaded region only and completely ignore the un- 
shadeff regions, Ki'pirc 1.2. The tnethod was used by Etkln and 
Woodward to obtain an api>rt)xiinate solution for a flat delta 
wing with subsonic leading edges. They also gave a procedure 
for arbitrary dowinvash disfcr ibutions . For the case of a jilane 
delta wing, it wjis shown by comparison with exact solutions, 
the error intuirrod by in-glccting the contributions after the 
third reflection, is negligible. 

For sharp subsonic leading edges, the attached flow 
theory predicts eitln-r infinite or zero load at the edges, 
iiecause of the teruJency of the flow to separate v;ith the 
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infinite type of load at the Icadins; edf’ce, it is more usefnl 

to consider load distributions of the latter type. There is 

evidence that win'.>;s with droopei! edpes sust.ain to a ‘greater 

f>n 

decree, tlie thenretica! value of suction than plane win";s . 
Thus, ca'’iher<?d win.'TS may he oxpecteil to b<! better than T>laiic 
winps from tiiis view point. in what folloivs, we discuss eases 
of attached flow and zero load nl on’; the leadtua, edi’es at 
thfi dcsi''ni conditions unless otherwise specified, 

1 <h>ni ca 1 i n»: Thco r y 

The pj'f)hleri of w'in-s of arbitrary planrorr.! with cott- 

r.letcly subsonic lead in, e; ed/^cs has generally not been solved. 

Tile flat d(?lta wing with subsonic edges has been solved app- 

2C) 

roxiinatcly by Kthin and P'oodward as inentionod earlier. 

llov/over, more elegant and closed form solutions for the delta 

wing have been obtained using conical flow theoryoriginany 

1 

<itio to hnsemann . Conical flo^v is <Jefincd ns a flow field 
in w’hich a'l the velocity components, and hence the pressure 
and slope of streamlines, arc constant along rays emanating 
from a point (the wing soex in this easel. Conical eamlicr is 
the name .given to that distribution of angle of attach which 
produces conical flow over a win;; surface. With this restri- 
ction, the three diraensional i’jandtl -Cl auert epur-tlon can 
Ite reduced to a two dimensional iaplace equation and the ivel 1 
devclojxMl method of complex variables may be used. 



Early invest! cantors usinc^ conical flo-w assuripti nns 

also made the 3 ov.' nsn(?ct ratio assumption thereby rendcrine; 

their solutions indeperulent of Mach nuriher and devoid ol 

the wave dra;r. . Both the flirect and the indirect liftin?^ 

problems wore studied, ProhaTsly the earliest of such atterrpts 

was that of loncs^* which reso] tod in his famous low aspect 

ratio theory and ,‘'ave a sinuile and ele'^ant analyticnl solu- 

1 C‘ 

tion, Thwier shiilar assunptions Sliaw in his studies of 
leadin.": edyc contr<»ls on a delta wins; "jnve the solution for 
the case wlieit the central portion of the v.-in,'j; is strui-dit 
and tlie outhoard ])orti<in drooped. The downwash distribution 
i!i his case was of the form 

't&'/U - a ^ ^ ^ 111 ^ V 

=. a ^ b c^-^h ^Ivl ^ ^ 

where n =; n, and w is the z-per turbation velocity, U the frec- 
stream vifloclty, the conical coordinate, a, b, n and ^ 

O 

are constants. I’he work xvas extended by Brebner"" for the 

9 

cases n = 1 and 2 and by Weber , who in ad<litlnn, considered 
the constant a, b and ^ n.s functions of the streariwise 
coorti i note , x. She was, therefore, able to extend the method 
to cases of arbitiary planforr; ;<nd downwash d i s tr il>uti.on for 
very slender ivitui^s. ’fowever, tliis kind of extension is tm’vli- 
c i t i n any low aspect ratio theory as .'yiven by Jones. I'urthej , 
the liaslc low aspect ratio S‘»lutioj\s arc conical , 
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'» 4 It! 

In the above work/"’ *”* the bnimiinry cemdition 

on the win”; snrl’ace w;is satisfied on an avixiliary plane, 

usually the silane z = O. fisinA the uhove assumptions but with 

the surface bosindary cf*ndition hcin" satisfied on the win" 

r, 

casitber surface, Smith solved for thr> case When the spanvvise 

3 

camber line is a circular arc. Cooke , in the sn^te spirit, 
extended it to tiie case wlicn the camber line has a striii'rht 
centra] ''ortion witli drooped edpres, a configuration obtained 
by start inrr f rnin a circle thron'Tih several traus fornuiti ons 
orlyina 1 1 y (Inc to Maskei!. Sjiith’s res ■! ts come out as a 


special cuiso . 

The direct inohieri of specify in'!; the loail d i stri hut Lou 
and dc tor'nininy the spanwise camber was solved by Smith 

(i 

and ;'an,':i;ler , uj^ain us in," the low iispect ratio assuiJi[»tion 
aa<l satisfy in*’; the surface hojmdary condition on “the plane 
z = 0. They could relate the twist distribution for different 
values of the slcnrierness parameter, but carryin" the same 
spaiiwiso distr ibution of load, and hence considerable atten-r 
tion was paid to the slender winp; solution. In particular, 
i t \?as shown that for slender wlns;s, the camber cou'l d be 
<lesi"ned to "ive lift dependent dra." factors as close to 
iinity - the tlu'oretieal minimum for wines sheddin." substan- 
tlaily niaiie vortex sheets - as <!ne vislicd. Thus, a win*?; 
canbered to have an attachment line alojc: th<> leadin" edye 
a t i ts (ic'S if'n Cj could he, ivxpx'Ctcd to show a useful reduction 



in lift depentlent drag compared with the sharp edged uucam- 
hcrcd wing, having at all C / <>, separation at the leading 

■Ls 

edge. 

A not very useful feature of Sintlt and Tangier’s work 
is the wavy cynber shapes that result duo to the representa- 
tion of the itotcntial function by an infinite series. As tJic 
niirobor of terms increases in the series, the lift dependent 
di'Og factor <Iccreasos aiid the can1«er lino contains more, hit 
less pronounced, vaves . A chorihviso wavcinoss in the ipreassuri 
li istribution givi's rise to the undesirahle feature of having 
several cor?nrcss i ons on the wing. 

The low asp<'ct rfstio assunotion appears tf! give ren- 
sonahlc; res 'Its when auplied to delta wings in supersonic 
flow. I’or cxan»»lc, it predicts the presstires near the apex 
of a plane delta wing with gooii accuracy anfi the *sl)apes of 
ihe spajiwisc pressure d i.s fcr i hn tion arc generany iji good 
(piaJitativo agreetient with exp.eri mental results, although 
errors in scale exist which inerfase as tiic trailing edge is 
approacliod ’ . however, the iow aspect rati<i assuriptlons do 
not dlstlngnish 1)etween suhsonic and supersonic sju-eds. In 
the suhsonic ran.go the thcciry fails to accoTUit for the I'utta 
condition at the tra i I i ng edge and causes an overcstimat ion 
of the ;)rcssurf's on the rear part of the vnng wi^ich gives 
rise to largo errors in the lift and pi tclving riojnent . 
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The low ospect rntio tlioory when a; piled to plane 
non-delta planXonns at incidenco shows that the lilt depends 
on the increase in width of the sections in a downstrean 

direction. Tlie sections behind the section of maxiinunt ^vidth 

develop no lift. i-Tirther, the spanwise loading is independent 
of tiie planform and <]ei>ends not on the area hiit on the width 

of the wing. A rectangular planforn tvould thtis have all the 

lift concentrated at the leading ed.go. 

The delta'wing prohlon, without the low aspect ratio 
assnini>tion, has been attenoted anong others by Vhildwin^'^, 

j— JWI. 

Tsien , Holla and in tlie present v;ork. The first two sitecify 
the pressure distribution as an infinite series of certain 
pressure functions and solve for th<‘ camber shape. The work 
of !fo!la and th<' present dissertation specify the twist 
distribution and solve for the pressure fiistributjon. 

Tsien’ s work sJiowed that when a finite nundier of terns 
in t.'sc solution are taken, nuniinun drag 

slmpo, has an infinite value at the centreline. This undesir- 
able feature may he avfnided by adjusting the appropriate 
ennst.'ints so that tfic terns eoniributing to this feature 
are zero. Uy doing so it was found that the incr<ui.He in 
^^b^^L^opt not surprising, us the 

centiai part of the v;ing shar)e acts 1 i ko a fin with very 
sma i 1 erfect on the ptossure c! istribution, Tsien carried 
out the analysis for both full Icjvding edge and no 1 eading 
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odf^c suction. For the suction case he Couini the miuiniiri 
Cp to he almost identical to the plant' at tin* saiMc 

except when tfie lead in, ed.'te lies close to the Mach cone 

ty 

when nf)ticeablc inprovenents of about in C_/C^ over the 
flat plattj with snctioii is obtained. The suction ease Is 
not of ."rent it-e'ortance since in practice it is tllfficnlt 
to realise its full va.!m? due to s<*paration. 


The no-snetion case shows appreciaJiic improvement 
ovei" that of the correspond i n" flat p'late. In fact, the mini- 
mnn is only slightly higher thtni t!ie p’inir.nm of the suction 
case. Tlic use of conical <^a^!hcr will almost entiroty overeome 
the unfavourable effect of tin? lo:js of suction forc(*. 


7 

The indirect ;irohlen was con.sifiered by ifolla wiierein 

•■) ' > •> /<’> 

ije used conical twist d istri hntifjns of tl\e tyoc ( 1 -a'" “ } ’ ' “ 

fur n=l to 5, whore a is th<* slentlerness parameter and 'Yj 
tlie conical c -o rd i luite , and obtained the corres [lond i u" pre- 
ssure (!i s tr j !)u t ions . lie also arriveuJ at si"vilar conclusions 
regard i 3 i‘r. minisiurf dray, c.s those of Tsien. it was found in 
this worlc ajsd in all the others mentioned earlier, that in 
tiic no-suetjon case j luiducticm in diay is obtaineci by- droop i ng 
tlic outboard portion oi tin* win," so that the pressures acting 

on I’jat part g i v<‘ rise to a tlirust conjtonent. The present 

IMi 

wori; considers twist <1 istri but ions of tlie type -rj' for posi- 
tive integrai values of n. Tiie pressure tii.str ibuti ons are 
obtuinefl in a s«iries solution in increasing pow'ers of the 
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slewlerncys paraneter, a. These rostil ts ai"c snhseqnently 
us<!d to obtain rtininHin dra,-^ v/i 115*3 under c-fjnality aiHi inequa- 
lity constraints. 

An obvious extension to conical ilows would be cases 
where the lu^locity coinuonents vary aecordin/* t«) sorse powoi- 
of X, y and z ulon;* each ray. Such flows arc called quasi- 
conical flows and its velocity potential is described by 
hojiuipp' neons fonctinns of tiie coordinates x,y ansi z and esiiitain 
the conical flow field as a s'cclal class of defr.ree zo-'ro . The 
quns ! -conical flows, because of lesser rs'strietion on the 
flow prosperties wo’ild natsirally '*ive rise to rcinirtun: dra*,': 
win ':s that are acrndynatnica I ly siorc efficient cotMpared to 
conical winqs. 'ows'vor , tlis* snathsnnatieal problc’,! is consider- 

j 1 IP 

ably nore co!;pl Lcated . Lance ’ ^ ns-nnsv others hiss sstilisesJ 
this fiictbod to ftivl tin- properties oJ’ ivinTS arrived at by 
usinr* the first few tei-us of a certain (puis i -coni cal series 
for the doivnwash distri hut i on . These i-<;solts, alojK* with the 
conical flow resi.'lts of this d i SHcrtati on , nay be us('d to 
find tlie v.ud'c pro’uMties by usln-r the lift canoe E 1 at ion r:ethod 
and also to obtain .uerofiynasiica I jy clfieient win's, 

Sonu. ('xper ipjonta 1 resul ts "n conical v.'in’s are ri'ported 
in reference and 03. In le i'ercnc 0 (>3 , iviups are desiqner! 
us i in* direbJicr ’ s twist distribution lor u delta win,'.', of 
3; (iC!*r<‘C sct::i vertex anple and a spamvise (>i!iptie thicAuess 

d Istr i liu t ion . The corri'spond i n,'* plane (’ei ta ivi u"; Is reported 



in reference t>2 . The results show the possibility of obtain- 
ing a theoretical zero load condition al on'' the loading edge 
(at a slightly hi'^hor incidence than theoretically predicted^. 
At inci donees below tlie designed value, separation on the 
lover surface is suppressed while the deve lopinent of a tur- 
bulent bound.'iry layer delays separation on the top stirface. 
Comparison Ix'tween plane and cambered wings shmv that both 
l»ave the same lift curve sloyie and n<rarly erpial jlrag at the 
d 0 s i .g n I nc i d e n c o . 

In another set of exper Lments^^'^ on a plane, a biconvex 
and a conical delta wing of semiapex angle 37.5 degree it 
was found that spanwise cambered wings with drooped loading 
edges acliievo higher leading edge suction than a plane wing, 
specially at high subsonic and transonic speeds. The suction 
peak at tlie leading edge oi a plane delta wing is re-iluccd 
by a region of finite suction over the ionv'ard fhcin.g urea 
by the IntrodTiction of spamviso camber. The Piodified pressure 
distribution helps flow attachment at the leading e<lges at 
higher incidences. 

f • 4 Drag Minimization 

tuie of the central nroblems in aerodynamics is the 
determination of miniriut' drag wing shapes for a given plan- 
form and given lift. For delta and d(’l ta typt* winsTs it was 

r,f) 

shown that twist j)lays a more important role tlian cambej’ 7 
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In this resisect conical camhcr !:ay be considered to ,<:;ivc 

iisefnl dra," reductions. It was affaln shown by Tsien^ that 

superposin" a r)arabolic camber on a conloally catjiborcd win" 

aave additional relief in dra;". Cohen'^'^ usiny; the quasi- 

11 1 ■' 

conical resoits of Lance ' ’ obtaincwl lower values of draq 
conij)ared to c<‘rtain conical slianes but sf»me of her win"s have 
n]vturned leading ed"cs which are ’irnne to flow separation in 
a viscous medium. It is also not certain lhat a!l low dra" 
v.'in"S obtained by us ini the linearized potential flow theory 
will show similar roduetions lu dra" In a roa! fluid. 

Wiiat is then required is that in the abs<'uce of an 
adequate viscous theory, t!u* j>otcntial flow theory be tjscd 
to assimilate realistic desiri.n criteria by introducing fjuanti- 
tatlve and qualitative information, some of which may have 
notiuni? to do with aerodynamics, into the minimum drag prohlcji. 
It nay, for instance, be desirable to have a eO!".bore(j centre 
section of the wing so tiiat , at the design , its incidence 
is smaller than for the uncamht'red centre section with th»‘ 
same slope at the trailing edge* of the wing and thus possibly 
decrease the hciglit of tlie undercarriage. It would certainty 
hi? desirable to avoid upturned leading edges as in Cohen’s 
design. It might also lie useful to have a innnotoaicnny di'C- 
reasing incidence fror: the win-' root to tihe tip for favouraliie 
stalling cliaractoristics. I’rnm a nanufacturi iig point of view 
it might he useful to hiive a sl-raiglit spanwise central nortidu 



Ife 

anf! a dirioront contour in the outboard ro/^ion. br the win?? 
may bo required to follow a certain given contour over some 
part of the wing. 

Clearly, these orofjlems will arise in practice and 
the present tiieoretical design methods are not adequate. 
Mowever, tiraler restricted conditions several irmiortant 
results in drag minimization are availai>le in the linearized 
tiieory. Thus, it is sliown that the ineJuced drag is a minimu! 
in stibsonic flow if the spanwiso distribution of lift is 
elliptic. The induced dra*?, even otherwise, is IndepenCent 
of the details of the ohnrdwise pressure d is tri luition. Siihse- 
((uently, Munk*^' showed that the dra-x <Iuc to lift is a minii.unn 
if the downwa.sh in the conihincd flow field is a constant. 

Tlie combined flow field is the field obtained 5)y sunerim^ios i ng 
tlie flow fields duo to the wing in forwaril and in reverse 
flows, /.gain, it was sliown by .Tones' ’ ‘ that the drag is a 
inininun for a given volume il the pressure gradient in tini 
freo-stream direction is a constant, an<l for a ‘..iveii iiiaximu!:! 
thiolciif^ss if tin- (ire.ssure is constant in regions forward and 
aft of the ruixinuri thickness. hno notices tliat these results 
help us to iflcntify minimum drag shapes. Ilov-ever, except 
under very sinqile c i rcurn.s tances , they <1 not. tell us v/haf. 
the miniKuim value of the drag is or liow to achieve it. i'he 
corxbined flow field resul ts are valid for both suhsonic and 
siipi'rsoiii c siMxeds. The utility of the cor hi nod flow Jield 
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is that it reduces a three dimensional subsonic flow problem 
to a two-dimensional one and, in general, simplifies the 
calculations in the three dliiKUisional supersonic case. 

A method which has considerable design value was given 
32 

by Graham using the concept of orthogonal loadings, defined 
as loadings whose mutual interference drag is i^.cro. Thus, 
along with the velocity components, pressures, etc., the 
drag Is also superposablc In the linearized theory. Graham 
fvirther showed that two orthogonal loadings, each carrying 
positive lift, when combined in any fa.shton to support a given 
total lift, will give lesser drag than either load su?. porting 
the same total lift. 


The most common procedure, by lar, in tiesigning mini 

mum drag wings has been to use the classical Lagrange multi- 

S 0 T G G 

plier method. Thus Tsien , Smith and Mangier , liolla , Grant 

34 

Cohen and many others have used it to design conical or 

other types of wings for a given lift constraint. The results 

of the combined flow field have generally not been of great 

help and we are aware of only one afjpllcation whore Ginzel 
35 

and Multhopp have used the constant downwash result to 
design xTings in supersonic flows. The orthogonal loading 
method was used by Graham and his cownrkers'^^ for some simple 
shapes. It has also been tised In this dissertation. 


The linearized theory provides some interesting 
theorems rel ating the forward and reverse flows, which are 



generally reiorred to ns the flow reversal^ thcorens . Thus, 

20 

Mimk showed that for a lift distrUaitlou in. subsonic, 

steady flow the drag is unchanged by reversal of flow direc- 
tion. Karm^n"^^ and Hayes"^^* sliowed this to be valid for the 
supersonic case also. Several more theorems relating the 
lift curve slope, etc. were derived by others and finally 
they were systematized by HrseJl and h'ard^^’ and 
wliich also account for the leading edge suction, if any. 

T!ie theorem is stated simply as 





(i.l) 


where p and OC are the pressure and angle of attack distri- 
butions and the subscripts F and H, refer to foiwvard and 
reverse flows. The integration is ovei’ the entire ivtng plan- 
fonn. Tills result is vaJid for both steady anduustoady flows, 
at subsonic and supersonic speeds, subject to the refp.irement 
that tfie iCutta corulition is satisfied at subsonic trailing 
edges whenever tiiey appear. 


lisin'!; tlie above nentionod rosui ts , 'genera 1 i zations of 
some existing theorems are [iresented in Chapiter 4, %vluch 
siiould he of assistance In designing low drag -wings. Also, 
the relatively unused hut nevcrtlieless powerful miitliods of 
mathematical programming (Ml’) have been ado])ted for designing 
n'inimiini drag wings under real istic constraints (which can 
often be stated os a set of inequality constraints). Chapter a 



is cor;!)lf‘teJy flevotcMl to tin s aspect . in Clta liter 5 the 
inetiioh of o r tho'^nna I !oa{ltn'';s is placed on a irsoro forna! 
footiny and its importance in rornntlati jh^ the mtnlmuin drse’; 
]>robicm as a ’'P probi eni is shown. The dii l'erent methofls are 
compitroci in Cliapter 7, 

I . 5 feadimr Vortex Phenomenon 

Tiic low aspect I'atio ulnne delta win;.*:, wliieh has a 
sma i i lift-curve sloiie> iinist operate at hieii angles of attacl: 
at low svK'cds to sustain a given lift. However, e.vi'ii at 
moderate i no i. do-nc es , the flow at the leading edg(‘s se])arate. 
The boundary layers on t!ie outer parts of the upper and 
lower surfaces fiow outvairds towanis the leading oflgc- and 
i.iiere coalese, forning irer- shear layers. If the Ifiading 
e(!ges are sharp, the prianary senjiratiou is fixful along tiiom. 
Ml a delta or a delta typo wing, a single shear layer forms 
at each ieraiing <Mige, and assutnes uinlor tiic influence of 
its own voticity, tin* forn of a spiral vortex wfiicii lies 
above tlu' wing. This spiral vortex f'xerts con.s i f’erab I e in- 
fluence on the iressure d t s tr I hu t ion . i'he onc-onin'; flew ueml 
the vortex sheet down and reattnclies it to t!u‘ wing, surface. 
The oncoming flow in the iidioard region of thi.s at taciu-unit 
line moves dii'cctly downstream, while tii*- liouiHlary layer 
which originates outboard of the attacluiient I i ne eneountor.s 
a ste<?i> pressure rise , as it passes outward.s utulor iiic vorti-x 
teavards the leiuifiig edge am! .separates tlius originatin/'-. a 




sccom! shciir !«'.vyer, I’i'curc" (’he so<i 1 e f this seconilury 

separiit ? on is sna! 1 juu! lias hriH! nc‘H ec ted i n ail thi'oro fc i n.'i i 
analysis out forward so far. niis typ<> of flow is found 
at. s’thsonic speens whcfi the sweepStacf' oJ tlii! ieadiir*' ed'^es 
is ahout o'; deyrens or rifire. A siniiiar s i. tuat i on occurs at 
Hiii-ersonie speeds when tfie leading (Hl.yes arc subsonic, the 
flow roattachinE witli inci easo in ^iach nurdier'* ’ . 


f'he leading cdgi? nhouoiuenon has attracted a yreat deal 
of attention iji recent yiaurs hiuiause only type of 

stable flow can he jsiai ata i net! over a wide ranpe of attitudes 
and "“ach nurthors ivitiiout the need for f i ow coiitroj (b’vi cu.s , 

The riiieiiortenon Is associated v;ith su iistantlal iiicreases in 
the !irt and tiie li rt-cur\u; slope, ’)o th of whlcl! tend to 
reduce t»ie need for hi‘r.h lift devices and hieh ineidenees 
fiurinE talce-off e.nd Inniliny. fh i s ?;;ode lof sepaiuition is 

it 

different fror’. the conventional houndary ittyer sie»aration 
in tl'.e souse t!iat in this case tlse incj’ease in drar. Is not 
oroi) 1 hi t i vel y lit 'h as in th(' case of h ountlary iryer sep;! ra I i o n 
Hie ,l callin'’; effye- vorti'X ■phenorieiion is thus an attractive 
Huhst. i tutf* for hiath lift devices. 


The phenor.cnoi! is undouhtedly due 
plays an i,ry>ortaut role in the hoaudury 
siirfacc, in the initial deve I opiueiit of 
and in a srsa L 1 sah-eore whi ch is found 
spiral vortex. Vise viscou.s s;th~core is 


to V i .scf;s i. ty v.ii i cii 
1 irvers cu tiu' w’ i n 
l.fie ih et- ;diea,f iayi'i s 
i n ih<- t-onl r*- of tf le 
ty ■ ^ i c-'!. * ! y no ;i:o re tlsuii 
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Sf-i O'J ill® ter o.f ‘feise woFtex*^’^ ®smI ®®.y h^^T^■ ax'i.fi-i ’S'e!®- 

citJfts tWG t-G tfjxe® il«a tije tT&m strmm mlHe, l',4srasbaw*'s'^’"* 
iiS’resfctfatS.#3s» stew ttet tb© total pressare In Use 
enr© tlw i?ot’t©x is os^siaider&biy l&^^v t-tmm tfe© s^jiTormtl- 
iitiS; streams ®-b<I ham its rslf^iisisr!; vain:© at the W'rtftx 
fiti® fast »ay 1« «*8ei to #!X|a€firi®eatally d<<=te2’:'ti'ri© tfse %-orl;-ex 
centre., ftatsixle’ ol these region® the flow f tel ft Sfipcar® t«. 
i3®p-aB:i® l:l'tt-le o® tbe I'eyfiftlSs aiiolj.-er aail t-fje flaw P6®y Iscf 
as®s»cx1: to fee In'-viscifto 

li'-Me to aeparatioRs t-tee pi'ess?/re peaJss ImMcaieS fej/ 

■the a-itsselisc'fsi ilfm theory' are smt dexclepe^s. hut- peafea are 
«feTel»|sed £i3!l»oa.r4 ®f ibe leadlisg «4ge®« i.xper iE»entaii|' it 
ba® lice» ©lseerw«l .tor fjelts tteat t-fee IxsiigltetfisMSl centre 

af tis© lescSlssg e4ge sssetloia,. tlic fwrtostle.l fl 0 ^ aeracfl^issB^ie 
cesitre'j, swi tbe aere’eljs^ie eeiitre witli leadlrii.g e«}gc- nortex; 
are s;»pro:s:l*sat®ly tfoe sme.. It ■af«?ild felKi® tlist tfee 

reiiist'r'ibfiii®©. tiie pstestlal tlOfS' part t-otB-i lift. i» 
prisarl ly aparwia.#. 

4T 

Mirfe-potricls*® ' esperteenta atewesS t-tet the path trsc-iKi 
■|^ tlie ■wnrtex: ecKtr«? :stew a nimfj a® ite S.aeii!e»e« is 
t® ieantJjf affeete-fl Isy tiie «?i&r4r»i sr?rle o# the 

ayjii feecases 'jwre ©?irw»S e»ft ® I #»©■*■ to the axto. ;a® t&c d-i- 
.fetrsfr®*! is I Bscressei", iiartfaer* fer » rtewfeie ers>s®-*»,«etl'r>fi t1'r.e 
wrtex centre is closer to the ilsa® for a fe-lcoBTex 

eres^s-sectfais lafttia®, toe satac leaisIS®?: esiee and 4!^4e4rsi an-r;l<;'-a'. 



Ai sri , tli(‘ force cx)fi'f ic'j <.‘3!t is ssnal ler for rhoj.'.bic 

cross— see tif>ria thoi for biconvex cr<»Hs— stjoti ous , and 
decreasi's as tiio dihf;dra! an.<xic inoreusc.es. 


^issoc i a ted wi th this is the pliomrieisosi of vort.ex 
Tji-eahilO’.vn f irst observed by feckha:-; and A tl; i ns-'a'* '.vhich can 
bo a iirritiji'': factor on the f> ft a tin-- aftitnr'es of slonf'er 
wMiyed it i rc ra f t . fh<! jioss i !> i 1 i ty fi f its occurrence li* the waive 
of a larr-e aircraft rsay lie fiazardnns in dense air traffic. 

■ freaV. dov.-n nay oecur earlier ovi'r oth* siiie of tiio ■win.'’’ .«-ivin'" 
rise to usioi t'asairt ro 1 M n'r a.nd nitcltinr’: vsntions and loss cs !' 
lift. '•’his ;^)hc rionennn is !>e i re’- videty st'idied ’noth exaerineii- 
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fio'tt’ dovais trefii ' appears to drive the fiow i!!i the r u ;>s ts-t;aiu , 
ndvile a redticti on of vsressure by suction will iir.evc' i f. further 
downstreno and (.'ventua I I y elij’inote it a ! 10 ." t' itn.n- . riie brer'-- 

(bnvn SOS' > .r i H i ny I y seens te depend little ftn the h'r^ynolfhv 
nmiiljcr i»o til for .la-’iirnir iUKi turbu'l en t %’ortex cores u -est rea:-; 
of hrcodvdown. 'fh«r ef fec t of <;or.’prf’;''S i b i 1 ■[ t.y see!’ ,s not to li;i\'e 
fieejs s t.ufi i ed i n eetail i. II t'uc ! i teratu r e . 



S«vcra.! [m (,<'nt;ial fJow inodcls have been cobh true ted 
tc) nredict tiieorotlcal ly tht' aerodyiiarsic cfairoctcristi.es ef 
tlie icatJin'r edeje vortex I'hcisosiienon * These may ho class! Tied 
lute the Tolltnvin" six cate-veri t's : 

i. Cross r 1 ow theory — this is 011 erpirical proctalurc duo 
4 2 

to dotz based on the assui“ption tljat the la torn t cross 
now of the winrr becomes iiceortnnt if the asia'Ct ratio is 
very snail. The nor!-l inear tern for tiie force coo i'fi.cients 
lesol ts nost'ly froju tlie firnr. cr'efficient of this alr'ost two 
dinensional cross f I <mv . Calc iations arc veneralLy based r.n 
the cross flow of a r(’'Cta}r>'ul ar jslate with a flrav, coeff i c i. eat 
of 2.0. 


2. Ho Ilay-tlcrs ten theory - -.'his type of riodel was first 

40 

j'rosmsed by Ho I lay for a very low aspect ra.tio la-etanyiil a r 

% 

Mt tlu‘ instance <jJ* vnn-jairfiUiii , anfl ' a ciiord 

wise dis tr ilsfutloa ,of 'cf>nstaijt ntreiij^th ho yoi i'.ic,os. 


The vortices 

a r (' 

R tro tched ac roys 

the sna?i aiifi 1 

oirv/‘ ii t 

si*’e edves. 

The 

:f vortio'ys are 

no ionyi/'r ni 

■'O iC! a, n 


shod rearwards at oi) ^/2. ^ vliero (X in the eh nrr i 

to the wine; sfirJ'aee. herstei^''''. Hiths.eipHOitiy ,o^^ed, this. fi>r a 
delta wine hy saiifl i vi ili iiy It .ii'ito I rt,*o Lo.',uy!i l.iir ihi ts ■ 

ai'u'l 1h’oc:h!y niibro-iit ^ a thtofey .rer each- ■ of' i^ier:ears. 

■flio .rioii>‘'l iTiorir 'hrhia\" hori- ,oi' 4"^ wl th d .is . oi c‘dl qrii. te we i I 

a nh i.s pr 1 .nei pal ly ffae to- . tli o' rioi-e.*pl a.tw^ r-'-aj s i r Mo’i t j on of ' ' 

V^>r t 



3. Coiiieai tl'Sf?ory ■*' Tfiese nodela use the low aspoot ratio 

assumption hy vtotkiwi In the Treff ta p3 e^f: nti^ the golutiosis 

eay therefore bo oonsiderer! cordcai* An sttompt is made to 

represent the physical vortex pattern in this Trefftz plane 

by placing a oonoentrated vertex st the o<sre and joining it 

by an approp^t'tate feeding line to the ieadtiig edgo» Thus 
48 

Legendre who pioneered this method, places only two force 
free pi:rint vortices without the feeding line and arrives at 
non^nnltpie solutions as pointed out by Adarr.s^'^, This nosiel 
obvloualy violated tMa vortleity laws, brown and ’'icIiaC'T’’* 
modified the solution by connecting the point ■vortex t^-- tlif? 
loading edge with n straight feeding line thereby satisfying 
the vortioity laws* Mangier anti Smith'’ use Instead a curved 
feeding line, Tn these modeis the strength and pfsattion of 
ttie core vortex and the feeding line In the Trefftsi plnsio arc 
dotor-"firi<Ml by the hmmdary enuditions that -tlio nnrriol vcio- 
ettieo on t!io wing surface tuiiiioh and tho lt?adi,ng edge vortex 
shr,:ta be stream surfaces ami sustain no pressure diffortmoe^ 
ilowovor, the apiproxlmatiaisa made show a pressure differcnfe 
across the feeding sheet, so long as the isolated vortex CEm- 
ilnuos to grow 'In gtrcns.tli in the streaniflsc direction, fience 
an irit-cgroi torn of boundary eondltion is aiJplicd, liftrioly, 
the »ct .force o<5t1ng on the vortex aystoi- must vonish* The 
roH'iitg show quaJi tnilvc egreement but not good quantl tatlve 
agreonent . 
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4. K'nohonarm-Sqnire theory - KuthoHaar;^’'^ and later Sqnire*"^^^ 
used a wodel wherein the vertical dtsnlnoer'cnt of the vortex 
core was neglected, and the vortioity was aasuned concentra- 
ted in a thin sheet in the 2»'5 plane, which is assneied to give 
a piecewise continuons downwash distribution. 'Ilieir choice 
for the downwash wos 



° * h( « ic 


- 0 ^ 2 . , 
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(1.2) 


where , oC^ and'’^ .^ are constants and 'fj the conical c oor-linatc. 
Three 0f>iu31 tions are needed to evaluate the constants - Two 
of then are, the load at the loading edge niust vanisJt, and the 
average downwash nust equal the incidence, I . e. 

^ -i- oi-si . * ( 1 . 3 ) 

For the third condition ’'(ticheinann assumed a linear variation 
of withiCJC, such that e i for oC = a and ~ u for ac = ‘^/2, 
and solved the problem using the low aspect ratio assiimptlori. 
Sqistre re’^oved the low aspect ratio assumption end equated the 
lift developed by Ills theory for a a 0 (sJender wing') to that 
of Mangier ami S-r-ith"'*" as the tfiird condition, lor thin wings 
the results are good, 

%Ci 

5 . 'itiC: suction analogy --A thnri^ introduced by Vt, i.hanufe*"'' 

If based on an intuitive H«®<vq»tion timt llie* normal forc^ on ; 



tlie suj'lace rofj’iired itiOittiuiti the tiov' iUxMit the 

vortex is the Ha-ne as tho leadin';; e<J";o snotion force rcoMi rod 
to riainta, in attached flow ah?>i»t the leadin''' odfc in the 
ootential flow case. Tn s-thsonic fl oav , thir snetron f/>rce 
nay he ca ? ns i n,"; any of the. tvcl 1 os tahi i shed ti ftiri": 

snrlaoc theoi'les. i'or fvlinw (’(elta, urroiv , and diaremd 'lan- 
forms, correlatiorj of vs. cC v.lth exneriinents over a rumr 
of s"hs'>nlc tn Hor,ersoiMC ’:ach nninhers is excoilent, 

r> . Mastic! siniJarity theory - .’i i !m i 1 ar i ty enrin Krters 
(diat relate i,he vitri.ex deveioonent with Ineidetiee a.hnve 
slender, conical r.'n";,s \ver<; fouml by i as'ic 1 1 . Ms in."; these 

the (i<;velon!;tent of the flow fieid afonnd each of a whole; class 
of slejider v;in,''s cuji he jirtslicted, orovidcfl tin.' devein 'Ment 
of tin; field around on<* raiuiher fM' the class is (;•> ter!, i ned . 

i.’he vaiidit>' of ' a.sketi’s theory Avas chcf-iKh! Ci.'-.iiinst 
cxi)er ir-ients ])y Ki rhoatr i cl; ' * t>n a number i'f slender, conical 
AVin";.s Avi ti'j fljiferent cross-section shin-es, learliu"; e«l'‘;e anyies, 
and d'liedrai anylef-; at sohsonic speeds, AYhich de' oiis trated 
to a luuiiarlcal)! e i!(‘";r<‘e t!ie validitA' of Maske'll'.s tIuH?ry. 

Tin' nositinn of the vortex <n)re (suitahly scaled) akis found 
to he indenen.dcnt oi' the cross-section shay.e ntid to 

deoend only e>a its i ead i n"; an";ie and ei .'■':ej!(! ra I j sed 

incidence piirameter Avlu'ch i .s a function of Mu* win" "loouii* try . 
the re.si'lts arc exoc'Ctofi fei he vafif! Ji t .s'foersotuc soeeds a 1 S" . 

Tn t’'f> tirird ch.'i'rtf;r ive ujuaucsc a rvcA" r.-otii'l Ariiieh 
uses the Pol linwins ’ Slit: ti ors anaJo'’i:y auul a 'mod i f i uid K’tohe'naaii- 
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Squire an|troacii am] moso tlic ;>robloin as a miniutui;! dru>i; 
prehlein. The metliod {vives jvressaro distributions v.'hich 
a.^ree iairly v<'c i J with exo(*i'iitif‘nfcul rcsol ts fox' a ]>ifine 
delta xrin,". 

ft way he noted that Carafoii has also attaokecf the 
prohlon in which the tfitc’-ncss of the win/>: is acconnted for 
‘I'he analysis is nnite tedinns and jen''thy- "'e hav<!, in 7 o<‘ 
with the r<’st of thr» 1 I tera to re* , nt'yff'oled tlie thichness ef 
and assiitt’c t!ia,l; siHdi f'rre (;(,,!5 are se;»er aoKoh | e . 

^ bone 1 ltd i n"' 1 e!*a rl< s 

■f'. su!.:::a ri 7.0 , the; a«>rodyna' ;i c; ;nitl ’:eoj;!etiic [u’oof'i' i, i <• 

of eor.iealiy cav hcreil v-'in-'S with .sahsoijlc li'adin”; (‘dyes 

riyin’' at s ihersonie ?tach nniiif»ers ar<.‘ .yiven in Ohu;>ter I!. 

Their of X-d<;s i,'j;u perforinance at .stiufj and lar.yc incideMce.s 

« 

t'.re .''iven in Cli. .3. I’iie (.1 eoretlcal i:,o(Jci. yri'oo.seil for 

tliC lar",e ineid'noe Ct;se is vciified ayaiiust exoer i luentai 
data for a ylanc dfdlta winy and is fo'inrl yt^nerally to he in 
^>;ood ayrticnient for low asoect riikios. 

■file dra,": iiiiiiiiizatioii t'!<‘or('r;s are rt'viewed and exten- 
d<‘d in bha;'ter 4 T'er Jiftin": Wfn";s. Ch.anters 5 to 7 disens! 
diffr'rent laethfuJs oT obtainin'' I'ini'ira-, dra" witi’s snh.i oc t(‘d 
to f'qnal 1 ty anc! fneioialtty constraints wif.li tin' liolp of 
o r tho.yonai] londinys anri the a ■o'l i ea 1; i on of i:ia Mnsra 1 1 ca 1 
proyra’ rii ny . 


•It 11" 

feci 



1)1 


CIU\i'Ti;K 2 

CV' ,\'i cjA i i Y CA> ; •• '.j Ji; ,■ ; ! i j . ■ja ; ;t y- ; y 


i ''‘’vojrnjn^ Eotu'itinns 


i’he txirtlal rt i if crotitial cn->ntion <'ovr-rivin'’; tVic 
f .lnr,' past a thin wi nr.' In sitporsoiil o llfnr i,s .'';ivon by 




*2Z 


/ o 




v.liort' <f> is the tx'rtnrhatlon velocity potontlaJ, ’! is Lhc 
r,r('<>~siicur! Yaoh tnimhcr and {x,y,y/ arc the cartesian cnor- 
(Hnates. ' *r i tten .In eyi indri cai polar ooomi inatos , we have 




XX 


ri 


OQ 


1 

r 


i 


■'.'hei'e 

y = r cos , z = r sine (h .'3 

m 

If the flow is conioai one i;:ay use 

x/jgr ~ X r ^ = y 

nnfl write {'A *2) as 

(X^-O (--r. 

where % < i eharac i.eri sos the exterior of tJu' Mach con<> <11, 
the apex and % > I tl»e i.ntcri»>r of tin.« ''ach cone. If one 
e!u>os<'S X =•= ah^ forX^ I, then (2.5) reduc<>.s ir. thn l aplace 
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equation valid inside the Ilooh cone. 

% * • ( 2 . 6 ) 

If one chooses = cos^ for 'X<l, then {2,5) reduces to 
the wave equation valid outside the Mach cone. 

^ • (2.T) 

Tfere >') corresponds to the Mach cone and ^ = the cone 

axis. -^5= 0 oorrosfmnds to the T'aeh cone and Tr/2 the plane 

X a i). "’c asautne the wlnq to lie in tho reqlon x ^0, so 

that f = 0 for X <0, Fif^nre 2.i. 

The refliion of imediate interest is the wln^ surface 

lyliif? inside the Moch cone. The utirface condition requires 

tl'iat v<.flooltle8 nornal to the surface Eiuat vaniBh. An 

add i tionnl condition is required to hi end the flovu-ji inside 

and oiitaido the Mneii cone. This :ls obtained fr»ii the fact 

that ’ tfjcr«; exists a,qreenent as to X- =£ i for all derivatives 

of ij:' » dcfltir-tl In the region inside or outside the \‘,ach cone, 

14 * 

provided that there is agreenent for the function itself 
Since disturhenees oannot propagate u,)stroas!'j of the ’’ach 
cone » it foJlowa that if the frec-strenn is iiniforn, the 
perttirhafdnn velocity ])otcnttal p will he ascro on tiie Mach 
cone. 

To ohtaU's oxprossiofia reiatliig the T^'i't'irhution 
velocity conponents u, v, and av, the Irrotntionallty conditlion 
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is used vdiich inejuis tliat ■ ’ 

X (l!i + y (Iv + z (hv = r (jSxdu + cosS dv + slii6 dw; 
be an exact dU fereutial . Was can occur miiy if this 
expression is identically eenal to zero, with i'. , v, w bein- 
functions of % ami h. hence, 

^ %iln + cos <4 dv + siiiM <lv; 3 ‘,^2,6. 

At.tdris sta,'’o it is convenient to use complex varia- 
1)1 ea. Vi’e , tliereforc, define 

^ , 5 =: 6 - i^ 

) = h(0,^ } + iu'0,^ ) 

V( ^ ) + lv(0, ^ (2.**) 

W( ^ ; =: ^ ^ ... lw(», ^ } 

where the functions U, V ami v) are analytic fiuictions of ^ , 
tl:e real parts of which ai-e tlK; per ttirbat i on \M'lociti<ws u, v 
ami w ifspee t i v<j] >• , 

Vi til further ' ;ani pulations hetweeu (2.8) a.nd (2.h)' the- 
■Cot lowin'!; cfiuri, ions result 

~ J0 dU = (jV/cos5 = dV^siu^ (2. lit) 

which are known os the eoi’ipat ihili ty cone! i t i ons . 

1 11 or i )ic i pi 0 , if th(' boundary conditions are surficient 
to do ter; ; ine (' I ther I' or V', tlie rei'.ainin'i components can he 



round by ti.rc-yratiu-^. Thu into-ral/ioj, constani, is iv,„n.| 
from the condition on the .Vach cone oxccot for an iiun,ina,;y 
constant u-i.ich is of no use tn The case v.hen T is aiv. n 

over sone re.,ions and V on others is nnre u'i tidcul t to solve 
n rather ronojcte Uat oj; soiotions u-ith dil.feront kinds oi 
bonneary cnn<i i t i. r»n.s ur<‘ iiste<i in Joints ami Cohen*'^. 


it is TeiH'raily convenient to use a conJ-ori:i;i;i riyiresei.t 
ation for studying this orobjen. dor the piuticolar case 
ol a liltin" delta win."; with sui)soiirc leading'' <'<l"-es and at 
zero yaw, fiarnfoli riei-iveri a rf' la ti.onsii i i)ei,v.-een tiie m and v 


(.!0!':"nnf‘ri tjj 

of Velocity 

oji t!»o n 

'in>: s'u 

r < 

■iac<’ 


ui^) ^ 

uc„ 

, .Ufa 

n$ f 

d-id-CV,) 

%■ 

tan S’ y i~yj^ 

w 

'^0 

dri. 

2 •d‘ 

a 

*■£ 


X in j 

1 ^ 

y/;?; 




i 

' - 

p-% ' 



and 

~%Tio) = 

c f d ; 

, ,x,. 2 

Icift 'i> 

+ 't. a 

If 

Kik) J 

^ did- in^\ 

t dr, 


( 2 . 11 ) 


t / ' dri 

to d‘^( ^ 


Ir 


‘I KLk) J 

fr dri .0 




(i- Ci Mriif): ' d(p^ Sj: 
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where ^>1 Is the conical coordinate y/s{x). s{x) Is the local 
semi-span, U the free-stream velocity, a s tang, 

k » / 1-B , and Is the seml-opex angle. K and K are elliptic 
Integrals of the first and second kind respectively. 


The net pressure coefficient or load at any point Is 
given hy 


I,( ^ . Cp - cj . 4u/ll ( 2 . 13 ) 

where C is the pressure coefficient, + and - refer to the 

r 

top and hottoiQ surfaces respectively. On substituting for 
u In (2*13) we hove 


i (:i ) 


4 G 


1'^ '7i 






/ — vn — 


X dvi- . 
dril 


(2.14) 


The first term on the right is siJigular and corresponds 
to a plane triangular wing. The second tom! is the load dis- 
tribution that vanishes at the leading edges. For a given 
wing and design Mach number, zero load at the loading edge is 
realised when C„ is zero. This will be so only at one Incidence 

I) 

which la given by (3.12). 

For arbitrary' w( ), difficulties arise because of 
the logarithmic singularity and the term /l-o^ 7 J when a « 1 
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in tha Howeror, closed form solutions' arc potsuiisig 

for many oneos o:f intorcet. 


Lona 


2*^ Tavist Diatribt)Et.ifv 

The dovfiKmsij dUtribotion w( '« ), ,,uich sbaU itei.co- 
fortb call tho twist dl stribiitloa, used here to solve 


is ta 


*«> it X 15 j 


ifhich 'to the author* » kmwledea has not been soivod 
a'horo, «t*b»tltutln<5 ^or w('»] ) and «:!spttndlng the torci 


3 »!!■ 


1*') 


• 1/2 


In a binowlal series 
( :t-»“ r(2)-l/2 . Jgn “ 


=" -1=“. a . 1 


where 


and 


T (an.~l) iJ . 

‘^2n " iSiij M • '^o ^ 


(211) •» « 2,4 »e,.*...( 2 n) 
( 2 n-i)ll« 1.3.G (2n-l) 

in IJq# (2.14), we obtain 


L Cr\ ) 




Jl'V^ 


I, S, r 'in cfi 

L ^ ^ -^.v fl 7 /” 

77 rM=a 


^ / T / - / 7 ? 


( 2 4 i 6 ) 



-^ 2 k - X 'I 7=- •7.-7 


1 n f-c.” rn I, ! li'^ I)y nartw and aln >1 i I’y i n "- v o find 

i-Ckti) 

zfl-^ / , ^ 

i 2 k == 77,7 X ” 


•r_ 

i'A 


. x.'h k 2(^-0 


< ~o 


{ V. . J 7 


If(>n<‘n I,. '.Id.; ray lunv l»f’ axiua-sHod as 
^ Co 


iil) 


^ tan 0 


€'\ CQ 


f?0 




/ShS //-7 


+ (a, l-iim-t-'n) 

2 71 


(2 .IH 


viic ('uiist.aiH. C ni i;:!.!::; V'liicn dot'.crtM IH'S Uk' Htranr'i.h of 
o 

!,Itf Ic.iiliii"’ H i 11".' I f fi f.y is Mvon !>y 

/* 

C f i 


tuin V / 

kn^ 

BCk) 

eLk) 

00 


sy f 

/ 7 Jr 

1... ■ 

-'o / - 

ns} 


n=# 


2n C7 


Of ft 

A Kik ) 2_ ^ ^ 

/f nrz't 




■f. . . Srj - 


X f (I- sin f) 

*^0.' 



f srt-l .{(i 

* ii-t ^ Jo '“ <i’P H 

in nnlvin- l.!,ii. n.-ml to solvo i„r oortoin intotruls 
whitiJl M l‘<’ < ■ Vil t n;i t C'Cl 1 )f'tow, 

t 


f in-t 




, L - } 


R. 


in 


i ■? 






h a\ 


drj 


< >'( nuM iJi'' zs srn-s wo havo 


R 


m 


/ .• j -) . i -ly 

i v4 ^ ^ (i - Ct sin &) ^ do 


i D ( iml 3 \ r / 1 . iw/ ^ ./ , 


( * 3( ! ' 


’.’horo ‘! i H (h»- cnii;i!(<(,o :t<>Ui rune lion end i’ (,lio hype r-^('(i!;;,- !: c • e 
i'une i.i oj! . 

Hiwilarly i'er l.lio i ji !,er;i-a i a.ipeJir i li'-; ijj i:.<‘ iiilrci 
Lor' iu,Ji«i»'' Lh.e wiiuorc hrac!<c(;H in ;:’.J9; v.e have 


r< Jn-1 f, y 

Jo ^ I (r~a'sini>) ' d<p dr/ 


1 

in 


i . 


( 2.21 


V‘h(' r r« 


i: ’"■ j '0 


'i. yi . 




5 2 


1 2r} + l z\ 

2 S 2 ■’ ^ J 


( 2 . 22 ) 
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and flniilly for the in the last teria we fl 

/ dtt-i ^4‘ 


lid 


i ■'! J, J 4 J, 


■^Vhere 




I KuXj 


'•i-n J 


I 


‘^'9 


J I- u"'/-/ J y ^ 

■I K {' i 'i t f „i , 

A * J / ^ \ ~"\f I 2 J i (■<■ J 


(2.23) 


Tie nee (2.10) now heeoiriee 

^inS 


C 

^ u 




] ( I I - K(aj j 
£ (A It' 


4;^ r 




4 £{1 ) k (aj - f ) / 

^ knit '^‘'1 


n 


J ; / -'« // +■ '• 

Off 


I ( i-n/. A.'^^ +■ ~ f 


< ■** ) 
t-ik) j j 


(2.25) 

Eqe. (2.18) and (2.25) to|»othcr dotemine the pressure dis- 
tribution at any ohordwlse station for a conically cambered 
delta wtUK ^ith subsonic leading edges and n twist distri- 
bution given by (2.15'^. 
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The spaiwlse dlstrlhution of ntt is f^iven by 


. f 

c [' I ^ 


'’/? 1 


'- [ LL'j;) . 

C, Q ■' '/ ^ 7^ ' 


(2.26) 


On aubetltutliiK for L{ in (2.26) wo have 






/ j 4 4 


. .. / t:^ 




f. :,.;j;z-- -rj J 2- ^ 

Ij^ f ?7.v| 




#.ay fl J 

'6w '”1 ' 2 c •>«+«} 


^ 1 : J 


L r ^r„/rr- <,fev - , 

Z- zrj 


he k ‘“ 


(f 


■n-~i 


whor® 


..o 








1 / ~‘- H 


^ /-p .f h 


(a.2t) 


yr f. ^ r %(.k-i) »-(fc-’«-f) 


J ^ \ -- ry 

ak (ahO ^ ^ .£. I'l'y/J!/ ' 




(2.28) 



For the surface shape at any chordwise station, wo 
start from 


;y Z 


} 

u 


(2.29) 


and the sE-ooordlnote distribution 


5R •* xz (y/xtan£) « xz{'>] ) 


(2.30) 


so that 


I 

J lA. 


y(r/) ... '7 ^ ...c &C*J ) 


(2.31) 


and the sponwlso surface slop© Is 




I Z: 

^auS '^7 


(2.32) 


Intcfijrotlng (2*31) we obtain 




M - •? / dv 


(3.33) 


whore K la an arbitrary constant and shows that the wing 
shapes are non-unique to the extent that a dihedral angle 
does not affect the aerodynamic characteristics In the 
linearized tjioory. In further analysis K is put equal to 
zero, l.o. no dihedral. 


on solving (2.32) for the twist distribution (2.15) 


?k 


we have 



43 


and 

z('>f ) « xz( ) 

gives the orcUnato of the wing at any point on the wing. 

Vhott t!je ordliMite is aensured with respect to the 
wing netUaii surface , we have 



(n.35) 


(2.3n) 


The lift awl drag coefficients are calculated using the 
surface values of the pressure and the downwash. The lift 
dependent drag coefficient in, therefore, 



/ / v^j) (' 7 ; 



where Is the pressure distribution due to the twist 

distribution w,,. *» 

00 1 * 
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Now define 

h ^ 




I t-ii CD =<>/ 

I f’?-' <7J "<7 


then 


AJ 


= «» z u 

C -0 


( C «3 <»0 


Z Z (/, .• 


The expression for Is 


or 


4 = i , 


^ > O 




TT 


a 


'm^Q 


I I 


Ziynti) '^7 


^ ?S>i 5 


T 


^ 2 'Wi ^ 


2 m 


-ff 




Z (jm^ i ) 4 - j 


mt'C 


" n I o-v‘)''" 

jj/'-'O ' £/•>/ 


( 2 . 38 ) 


( 2 . 39 ) 


4 /qh^ ^ 


TT 


A ■ 2 ^^' 




2 (wi-L) ^ f 


Z. Z; 

j 7 „ V 


( 3 - 40 ) 
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whore 



(Zkf2)!f 4 


The expression for which 
onsc is 


^2.4i) 

corresponds to the Hat plate 


L 


ti Tan^ 

~~F(k) 


The expression for 



is 







(2.42) 


“ 0 . ,.. 7 - 

y ^ I k,«^o 

(2.43) 

Notice that 1. « -d, . 

1 lo 

The lift coefficient is it^iven by 


= J LCD dy} 

‘-t '^0 ' 

tXi 

L II 

i - o 

(2.44) 

Where 1^ are given by (2.40) and (2.42). 

The oolloctlon of reat lts (2.18), (2.25) , (2.27) , (2.35) , 
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{2,36t, (2,3lO and (2«44) gli't* the 
geoiactric characteristics. 


wecesgary aerodynamic arwl 


of Reaitltia 

In section (2.2), the oemdynamic and geometric pro- 
perties of a class of conically car^bered delta wings with 
snhsonlc leading edges were presented. The geoBctrlcal pro- 
perties are especially easy to visnalise when the wing has 
a single twist dietriTmtion of the type Thus, 

one may immediately sec that for n > o the wing will have 
drooped leading edges and with Increasing values of n. the 
central portion of the wing will become more and more flat. 

The droop becomes noticeable towards the outboard side and 
rapidly increases us it nears the loading edge. If one defint-e 
an arbitrary point “7^^ as the location whore z( - 7 .(n)« jf.j , 

and ^ is some specified constant, then for tx givesi K « 

■ * /o 

will move outboard rapidly with increasing ru 


The spanwise surface slope will be low in the central 
portion of the xiing but will rapidly increase as It progresses 
near the leading edge. A similar trend will be shown by, the 
angle of attach distribution, which wil,i be positive in the 
central portion and will move to negative values near the 
leading edge. There will he a continuous decrease in the 
angle of attack from the root to the tip, a property that 
alleviates tip stall. The droop near the loading edge will 
help the flow to turn amothly, and hence nay be expected lo 
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predict properties closer to the real situation. The' wings 
are also easy to manufacture and do not possess Iclnlcs. 

The aerodynanlc propertit?8 are a little more difficult 
to visual ijBe. V/c restrict the following dlsousfion to casos 
whore the load vanishes at the loading edge, if referenoo 
is mode to Tables 2.1 and 2,2 and Figuros 2,2a«*2,3d, it Is 
notieod that the wing proportios depend strongly on the 
anount of spimwise curvattiro for n given lift. Thus, with 
dooronslng onrvatnre (inoreaslng n) , tlio wing Inoldonoo 
doercaios and moves townrds the flat pinto value, and the 
suotion penH moves outboard and Inoroasos In magnUudo, Tho 
spanv.tse load distribution tends towards tho ideal oiH’ fcie 
distributton and progrosoivoly dooroasos the drog ooofflclont 
and the efflcleney factor and inereasos the h/D ratio* 

Tho spanwlse load distribution i iinlil«e the flat plate 
which shows infinite slope at the loading edge, has ssero 

I 

slope; the magnl tilde of the load distribution increases at 
tho centre I and decreases outboard near the loading edge as 
onmporod to the plane delta wing with suction. The effect 
is progressively less pronounced with increasing n* This 
kind of hehovioitr is to bo expected since the ongle of attack 

docrensos from the root to the tip and the wing tip Is le-ss 
loaded compared to a plane delta wing. 'Hils gives bettor tip 
otall ohorocterlstlcs. 

The suction peak which inoreosos In magnitude and moves 
outboard with Increasing n, nets on the foirrnrd facing area 
of the drooped edge and thereby produces a thrust force which 
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TAIIIJ': 2.1 

; OVFFAiJ. 

AFFODYNAriC CirAliACTiJElSTICS 


• 

or iwsic 

SHAPES 


Tin" 

(n) 

s 

l/d 

e . 

a 

i 

.80109 

91 .05 

1.371 

.032 0.10 


.00008 

101,52 

1 .238 

.024 

3 

,00004 

100.40 

1,180 

.021 

4 

.00091 

109.4 

1.149 

.020 

Flnt Plat© 

.OUiU'SS 

121 .18 

1 .034 

.010 

(t7lth Suction) 

Flat Pint© 

(tt/o Suction) 

.00102 

61 .70 

2.038 

.016 

1 

.00151 

00.31 

1.895 

.035 

'i 

.DO 140 

71.38 

1 .700 

.027 

3 

.00135 

73.81 

1.702 

.025 

4 

.00132 

75.24 

1.070 

.023 

Flat Plot© 

.00124 

80.76 

1.555 

.019 

(with Suction) 

Flat Plate 

^ ^ V 

.00193 

51.82 

2,421 

.01$ 


(v//o Suction) 


1 

.00182 

54.89 

2.203 

.037 «•: .70 

n 

.(;0172 

58.07 

2.164 

.')29 

3 

.00168 

59.57 

2.109 

.027 

4 

.00165 

60.44 

2.079 

.025 

Flat Plate 
fr.'i.tli Suction'' 

.00157 

63.56 

1 .974 

,021 

Flat Plato 
(w/o Suction) 

.00214 

46.70 

2.602 

.021 

c s efficiency 

factor 

= 'S. "‘o 

= V’ln," Irsci flora 


tanS - l.f 
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amp geometrio firrA»Ar!>ygpf«mYr.ci 

9F iSASIC £'l\?E£ ”'™* 


® 0 « 50 


wing 

in ) 

1 

id ) 

Cl 



zd ) tmB 

/ 3 2 \ (j 

i J 

i 

.00 

.842 

.842 

.350 

-.000 

-.024 


0.20 

.891 

,794 

.327 

-.024 

-,238 


0.40 

i.oie 

.653 

.266 

-.095 

-.475 


0.60 

1.104 

.438 

. 13 ? 

-.214 

-.713 


0.70 

1.203 

.313 

.059 

-.291 

-.832 


£.80 

1.174 

.187 

-.030 

-.380 

-.961 


0.90 

.992 

.072 

-.131 

-.481 

- 1 .070 


(‘.96 

.766 

.027 

-.186 

-.636 

- 1.129 


1 .00 

.000 

.000 

-.239 

-.689 

- 1.185 

2 

0 .00 

. .700 

.760 

..273 

-.000 

-.000 


0,20 

.788 

.736 

.272 

-^000 

-.006 


0.40 

.891 

.660 

,.308 

-.005 

-.060 


0.60 

1.109 

.487 

.197 

-.025 

-.168 


0.70 

1.248 

.372 

. .133 

-.047 

-.267 


0.80 

1 .366 

.239 

.034 

-.080 

-.399 


0.90 

1.206 

.100 

-.110 

-.128 

-.508 


0 .95 

1.068 

.038 

-.203 

-.159 

-.607 


1.00 

0,000 

.000 

-.300 

-.103 

-.768 




mi 


TABLK 2.2 ( Contrt .) 
a = 0.50 




4 


1 

Lj^J 

Q 


iX /tjHF 

z 

Ci, 

/ sz \ ianS' 

> 7 - 

'L 

■n .00 

jnti 

,725 

.247 

-.Ono 

-.000 

.20 

.747 

.705 

.247 

~ .(*{(■:> 

-.000 

f ,40 

.827 

.638 

.244 

-.■’Ol 

-.008 

'. f )0 

1 .035 

. 5».»0 

.218 

~.fH 6 

-.057 

0 . 7(5 . 

1 ,215 

.405 

.174 

■ -.014 

-.124 

(.'♦BO 

1.425 

.275 

.085 



0.90 

1,514 

i <^1 n 

-.081 

-.065 

-.437 

0 .95 

1 . 32 B 

.049 

~. 2(>6 

-.091 

-.672 

,l .00 

.!lOO 

.000 

-.352 

-.121 

-.722 


0 .00 

.706 

.706 

.233 

-.000 

-.000 

0.20 

.725 

,688 

.233 

-.000 


0 , 4(1 

.794 

,629 

.233 

- .*000 

-.001 

t .60 

.978 

.514 

oop 

a 1 , £.f 

-.002 

“-.021 

0 ,70 

1 .164 

.423 

.196 

-.005 

— , (}62 

0 .80 

1.435 

. 3 ( H > 

.123 

-.016 

-.157 

0 .90 

1 .066 

.142 

-.049 

-.040 

-.368 

(• .95 

1 .550 

.058 

— « 20 1 

-.062 

-.523 

!. . 0 (! 

,000 

.000 

-. 3 '^7 

-.091 

-.724 
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heirs to reduce drag, lor 0=4, the Is substontlally close 

to th© ilnli pliit© veIu© with ©uotlon. For © ti it 

approaohoB the flat pint© dra?? faster with tnerenslng valaos 
of a» The leadlMiii; ©di^c thrnst obtainable in a plane delta 
wing da© to leading edge auction, ie obtained in the cae© of a 
oaiiihored wing with a leading edge droop. Thus, the caetbered 
wtngB, when compared with a flat delta neglecting leading edge 
thruBt, are definitely esiperlor aerodyneniically« 

Conically cambered wings may show some advantages 
over other shapoo under off-doslgn conditions. Thus, at any 
supersonic speed, the lofid will vanish at the lending edge 
at BOW© incidence and, therefore, the attachment line is 
unlikely to lie on the upper surface for part of the wing and 
on the lower siirfaco for the remainder, even away from the 
design Incidence and Mach number. 

Conical camber may be restrictive, since a control 
over the pitching noment Is not possible, 'fh© cliordwisc 
dis tribution of load la 1 ineur and tho acrodynfimic centre is 
always at the centre of area. This may be rcfncdied by iislrw; 
a parabolic emtber or some of the results of l.anoo . 

However, the simplicity of the surface shepes may be lost. 

l''oi* wings lying near the axis of the Mach cone, the 
wave drag duo to lift 1« a relatively small part of the 
dcpentlcnt drag end conical camber may he fruitfully used. 



S4 

For wins® with leading edge® close to the Mach cone, when 
the wav© drag due to lift exceeds the vortex drag, conical 
caraher will not be very useful. 

It 1® possible by choosing higher powers of -yj that one 
may ooroe a® close to the flat plate value as desired, however, 
the eparKwisc surface slopes and the magnitude of the tip 
incidenoc will become very large and the small perturbation 
aBSumption® arc likely to be violated. 

From a practical point of view it would be interesting 
to find upto what value of n the wings would be useful. 

Invisoid potential flow theory provides no answer. The matter, 
in practice, 1® closely cormccted to the design C, , 11 and the 
various basic shape® that are combined to produce the wing 
since they determine the adverse pressure gradient® which the 
boundary layer must overcome. They will also determine, to a 
certain extent, the state of the boundary layer. 


From the analytical results mentioned in Section 2*2, 


the low aspect ratio result® may be obtained by putting aaO. 
These result® may then be used, as done by Weber, to design 
wing* of arbitrary planforms and downwaah distribution by takir 


the constants a® functions 


of the streamwtae coordinate, x. 
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in the integral remains the same. We can, therefore, for a 
given pressure rtistrlbntlon, find the downwash distribution 
reqtiired ot n given Mach number and apex angle. This is 
easy Iti our case since w Is expressed in even powers of ^ and 
1 / expanded in similar lowers for a < i, and 
the results of section 2.2 may be applied. An interesting 
feature is noticed from Table 2.3 that for a given value of n 
and C, the ©panwiae pressure distribution and the 8panwis«j 
lift distribution differ little for different values of a, 
although the wing shapes and twist distributions are quite 
different# Thus one would conclude that the wing root bending 
moment would be about the sane for these cases. 
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TABLE 2,3 S COMPART SfiN (■)' WINGS AT n i FFl-:RENT VAI.JTflS 

OF ’a* FOR ttr: msic shapes* 


n s=, 


n 




l-CT)/ 





1 

a~oi 

A = o-S 

CL=o7 

a=oi 

a=5'‘^ 

a~o-i 

0,00 

.T64 

.760 

.757 

.764 

.760 

.757 

O.20 

,791 

#788 

.783 

,738 

,735 

.731 

0.40 

.807 . 

,891 

.884 

.651 

.650 

, 648 

(V.60 

1 .115 

1.109 

1 .099 

,485 

.487 

.489 

0 .TO 

1.251 

1.248 

1.242 

.369 

.372 

.376 

0.80 

1.350 

1,356 

1,363 

.236 

.239 

.243 

0,90 

1.276 

1.296 

1 .324 

.098 

.100 

.103 

0.95 

1 .045 

1.068 

1.104 

.037 

,038 

,040 

1,00 

.000 

.000 

.<H’0 

.000 

.000 

.000 





Z i &»S / x Cl 


1 



dr- 0-7 

(2=0-1 

a= 0 

a -0-7 


0 .00 

.241 

,273 

0 .20 

.240 

t}fn 

0,40 

.225 

,258 

( .60 

.159 

,197 

' f .70 

.089 

.133 

0.80 

•» .019 

.034 

0,90 

-.175 

,110 

0.95 

-.276 

-.203 

1 .OO 

-.382 

-.300 


.293 

c 

o 

1 

-.0t)0 

-,00() 

.293 

-.000 

-.000 

-.000 

.280 

-.005 

-.006 

-.004 

.226 

-.027 

-.025 

-.023 

.1C57 

-,05i 

-.047 

-.042 

.078 

-.087 

-.080 

-.07::.' 

.052 

-.139 

-.128 

-,li5 

.136 

-.172 

-.150 

-.143 

.223 

-.208 

-.193 

-,.173 


^ I A 4 wi rt iM. Blow* wli©rdis results of 

Tslew^^nolla 0to!! 0how murked changes. Ilenoe, ofl’-tJesign 
Tslen, lloJia, exo. »mm to fee rauoli better 

«ax Ho.. 





TABLO 2.3 

n « 

( Contd . ) 

4 
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1 


^-( 1 ) /q 





fl “■ O'l 

a-“o*9 

a = ot 

^ 0‘l 


a=o‘7 

0.00 

.707 

.706 

.704 

.707 

.706 

.704 

0,20 

.720 

.725 

.723 

.689 

.688 

.086 

0,40 

.796 

.794 

.791 

.030 

.629 

.628 

0.60 

,982 

,978 

.971 

.514 

.514 

,515 

O.TO 

i .169 

1.164 

1.156 

.422 

.423 

.425 

0,80 

1,437 

1 .435 

1.430 

.298 

. mo 

.303 

0,00 

1,657 

1.668 

1.685 

.140 

.142 

.145 

0.95 

1,529 

1.550 

1.584 

.057 

.058 

.060 

1.00 

.000 

.000 

.000 

,000 

.000 

,000 






; ’^'>1 S 1 X. C-i 



£l o ' 1 


(X - 0 * 7 

a~o-i 

C'-S 

« ~ 0*7 

0 ,00 

.202 

.233 

,255 

-.000 

-.000 

-.000 

0.20 

.202 

,253 

.255 

-.000 

-.000 

-,000 

0.40 

.201 

.233 

.255 

-.000 

-.000 

-.000 

0.60 

.189 

rjjri 
% IM 

.245 

-.002 

-.002 

-.001 

0.70 

.160 

.196 

.222 

-.006 

- .005 

-.005 

0,80 

.080 

.123 

.158 

-.017 

-.016 

-.014 

0,90 

•* . 1 1 Oi 

-.049 

.007 

-.045 

-.040 

- .035 

0 ,95 

-.27 9 

-.201 

-.127 

-.069 

-.062 

-.055 

1 .00 

-,496 

-.397 

-.299 

-.100 

-. 0:91 

-.080 
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CHAPTER 3 

OFF DESIGN PFRPOrafANCE 

3.1 Ofi>-Deslgr> Calculattona ( am) 

th© off-fleslgn performance of a wing must be satis- 
factory to make it acceptable for use. About the flesign 
Inoldence anti Mach nuinbcr» it is, therefore, desirable to 
know how much penalty one pays if the incidence or the Mach 
number or both .change by small values. Of principal interest 
is the pressure distribution, the spunwlse loading and the 
lift and drag coefficients since these present a comprehensive 
picture of the wing’s aerodynamic efficiency. Linearised 
theory is assumed to be valid. 

3.1.1 Increase in angle of attack, 

The effect of an increment in Incidence ao <, can be 
obtained by stipcrlmposlng the pressure distribution of a 
flat delta at an incidence &o(.. The flat plate results are 
known, and the off-design results are given below. Hence 


^ ^ 'ii 


, , 4 /SuVl S 


{3.1) 


-■= CfCl) 


tCk) 


r,ci) n-7}.^ 


( 3 . 2 ) 



59 




2 ^ tanS 


and tho dra;; coefficient ie 

a .-- f' (/, + 1 r <K 

f •'“ ^ e[k) f,7^, ' 


dpj Ao<r J) d-Tj 


(3,3) 


- n 




C 4 -^ 

dfs 


jrr fart S 

£Ck) 


C^oi 


0(3 

> / Y " f: 


(in- 1 ) II 


L ^=i'H •7-'"-' ~ - e I 

wif 2'n (zr})!! -' 


(3.4) 


Whore the eubeorlpt »off» refers to off-design conditions 
and ’dcB* to the deelgn condition. 

I; 

3.1.2 Increase in ^faoh number, fiM 

A change In the Mach number implies a change in the 
slenderness parameter, a. Let us define 


(a * AO j ~ I (M+ mf. i ] 

The expression for the pressure dtstrlbutlori is 




(4 * J J0 ) 
^(Jn b //- 7'' 


f i<r 

/■■■■'■-•: r 7 


t a. f'*?! 


(3.5) 
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EaEpandtng In a lilnoi!iial aeries the term 

[t- (atoa) J iaV^- .... (3.7) 


where A a la the change in a due to AM change in the Maoh 
mmher and AC^ la the corresponding change in C^. 

on stihatltuting (3.7) into (3.G) and retaining terms 
containing upto ca" we obtain 




k S 

V ' 



.j 

' ^ 

7 

L"Kf 


■xn*f2. 

( a o a. 




■ iTf - 



j ^ ^ ^ 

*'■ ‘"ties r™' 

■ (a>fS J f-'f- 

^ (<2A<a 7- ^ 

TT 

' V r 

J 



(3,8) 

t i'l) - f'^i'UdeS ^ 1' C 

'n V '7 

L : ' krfS 


{>4 

X 2l -<• 

>f -- 1 

i ^ ^ J 1 

(33) 


and the lift coefficient may immMtately be written down as 



J 


t/0 


/-C7J„ <*1 

t 
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T + 7 6 -l.^z 

frtwh •<!- 


an ■>■ 3 


2^ 2.^ C,2 (>»t l-k ) 


Tho cxpresBion for AC^ is obtained as follows: 


(3,10) 


c. ^ ^ a 


l9±l 

elk) 


f- fS 

l-k- z 


za> K.[k) 


2_ ^2r> 


• ^ . 

-f 2 K(i.) Z < e(a') - s;., i / 2 , 

tr n-1 ( { 


/ , -v 

1 I K(k) - eik) 7 

^ I i »=i 


KLo!) - J, 


where 


» 0 4* ^ H 


(3.11) 


Now writing the elliptic Integrals K and H and the functions 
Ilf S and T in powers of a and retaining the first few terms 


in aa, we obtain after some algebra, 

«1 

- J. J - Z i) 


i 3 


£ i(k} ^ ”^4--n 



02 


r ^ (a f 

L 'If ^ 


2. 





( 3 . 12 ) 


where the functions B, K and E are defined in Chapter 2. 
Notice that the change in due to am is of the order of 
(Aa)^. 

The drag coefficient is 



4* 


r. ■ 

•»=/ 


(adfS + 
2W-#' 3 


^ 1.. '^k [ ' 7 -- (iitj~k*f) 2i 

ic:0 r' ^ ^ 




I {yi-k-i-i) 


J. 


( 3 . 13 ) 


3.2 Off-«l»eetgn at High Tn<tidencee 

At high incidences » leading edge vortices appear whioli 
modify the proeaure dletrlhutlon over the wing oonsideraWy 
and there la evidence that for a delta wing with sharp leading 
edges, the vortex development is almost conical In nature, 

Here following rCucheraann*^'* and Squire we assume, that the 
leading edge vortices modify the flow In such a way thot tio; 
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fidditioiMil overfige downwash Is equal ■to the increftse In 
incidence from the design condition of zero load at the 
leading edge. Additional assumptions are that for this case 
the load at the leading edge remains zero and the additional 
wrtex lift Is equal In magnitude to the leading edge force 
os predicted hy the attached flow theory, acting perpendicular 
to the leading edge In the plane of the wing* This was shown 
to bo substantially true for flat delta and delta type wings 
by I’clhamus and Is assijined to be true here. Apart from 
these, the additional flow field due to the vortex system is 
supposed to ad.just itself so that It presents least drag for 
Itself. To solve the problem use Is mode of the conical 
solution obtained earlier In Chapter 2. 



The zero load condition is obtained by putting « 0 
in Eq. (2.26) 

- X’ I ( ^ ^ ^ ^ 

■ , /w : / C 
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rr 


'•L ) x-M - ] 


- O 


OT 


<s^ 

r ,.. T 6 A ^ 

4*0 Z. ^in 

r;a| 


(3.15) 


whore tho definition of Ag^ is ohTloue, 


The aiietlon analogy of Polharaus yields 


I 

4 : . I 




Af’f 


n 


SiV; 6 


C 


L 




fi-fct 


r/ftyjS (liay 


(f- 


'/j 


/ 2 r fertS £io( 


Sin 


e(k) 


or 


2 . ^u-» 


n~ 1 


r 5ec 5 7^ 

£•■( 1 ) X £~ik) 


a-^ 


4 - £ 


'iTivh "I 


(3.16) 


where C. and C. are the lift coefficients due to the 

‘ir 

leading edge vortex and that due to the attached potential 
flow theory for aa, respectively. 

The cocplct© problem would now read, if a finite 
number (N) of conical solutions are used, 


'•Inlmlzc (I ^ I 1 ^ 

t*r j*’' 


(3. if) 
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subjcot to 

N / 

{-o X 4 m / ) ■- - Zi£.t 




14 

2.» 

•n-f 




"in 


0 


(3.18) 


(3.10) 






Z'ri 




zr zio^ j- sec.'i) I i~ 

Vfti ^5(1)" 



(3.20) 


wlilch may be aolired for the unlcnoxm oonetants using the 

Lagrange multiplier method. The additional pressure 
distribution AL duo to the loading edge vortex may be 
obtained by introducing these constants in Kq. (2.18). 


nie total pressure distribution is 


and the total lift la 

4 - Q i ^ ^ 4 * 

‘■■Jfs ^att 

3.3 Mscusslon of llcsnlts 

For snail deviations In oz and M, the wings show good 
charootcrlstlcs as indicated by Table 3.1. Front a design 
condition of a » 0.5, C, * 0.10, and tan^ « 1, a 20^fc change 

M ■ 

in the angle of attack did not substantially alter the L/b 
ratio or the efficiency factor compared to a wing designed 
for the off-dcaign cases. Similarly for a lo^^ change in «, 


(3.21) 


(3.22) 



TABLn 3,1 OFF DESIGN PERFORJIAMJE OF CYNICALLY CAMBERl® DELTA 
IVINOS FOR n = 2 AMD 4 <s=0,5, tanS * l.o) 







tho wmo trend nne notlood. T:,e wlnss trny, therefore, be 

i ..Hf»lrS<'r«HJ J'or practical wae* 

I'or lar,"o deviatlona in a where leading edge vortices 

or© obe®rvcd+^, cxpOTlraental results are available only for 

t8 

a plan© delta hence theoretical results are coiii|jared with 
fcfiijsc* A reference to ilgure 3«1 shows good agreeraent), in 
genoraJ , with the experlnsentnl pressure dletributlonB for low 
aapoct ratio©* The comparison for the cose is”, 14 

shows mnrhcd deviations. This is the largest incidence shown 
and Indtcfttos that at such high Inoldenco, the n«n-lltiear 
ttiiohncHS effect Is probably very Important, a fact neglected 
in nur msalysls. It is also possible that secondary vortices 
have developed and have influenced the preseur© distribution. 
This fact Is also not accounted for in the arsalysls. 


* Liir/*'© iliivlatlon® are used in the sense that leading edge 
vortloeH have devoiopod and have a significant effect. 
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SOME TIIEORIJMS CONCERNING DRAG MINIMIiiATION 
4-1 iBtroduotion 

One of the oentrel problems in aerodynamics has been 
the finding of minimum drag oonflguratlons for wings, bodies 
and tholr oorablnatlons. The literature in the past fifty 
yearn hns evolved theorems that generally identify a lalnlmum 
drag oonflgwratjlon, but seldom say how to orrive at it. In 
deriving the nlnlmiiw drag theorems the combined flow field 
concept and the reverse flow theorems hove been widely used. 
The oowhtnod flow field | first Introduced hy Munt, Is defined 
aa the flow field obtained hy superimposing the disturbance 
velocities In forward imd reverse flow over n body hovlng a 
given dietrlbutlon of lift, sldeforce, end thlciaiess. The 
reverse flow theorem for wings enn be stated in a comrioct way 

wliero and are respectively the i>roBsurc and angle of 
attack distribution, and F and U refer to the forward end 
reverse flows, Tlie integration is over the wing surface. 

The tlieoioin is valid for subsonic, supersonic, steady and 
ttnsteady flow slttiations. 





TO 


III the followin'^ section we derive sone thcorenm ivhioh 
arc essentially f^enernllasatlons of those given by Rodriguez, 
et, ni, 'itie iiitoroatlng feature of these theorems Is that 
they help us in deciding how to proceed in order to reduce 
drag* The results derived are valid only if there is no 
leading edge suction, o,g, the eonloal wings of this disser- 
tation* The following definitions hold* 

4 2 ,2M. 

Thforeiti i. Any nomhlnatlon of t\fo loadings o(^ ) awl ) 

supliorttng 1) given to id lift will have less drag 
than either loiiding supporting the sane total lift 
provided 

' ida ! di ) ^'/ 4 . • 

\ 4 * a ^ 

Proof I Let ( ^ he two nonr^orthogoruil load- 

ings# Then after Graham, a pair of orthogonal loadings Is 
given hy 

pM B ( p4 ,^1 4 "" ^ 4 ^ 4 

K = 

where 

c « - 2 dfi I f<^/z ) 


* i («. 3 ) 


( 4 # 4 ) 



71 


Since J and «re orthogonal londin-s, they my be coabined 
in any nanner to reduce drag, provided each is supporting a 
positive llio, i.e, 

f k ciA >0 J ® 

■ (4.5) 

Uenoe it followa that 


- zdii j {du J *= c > - J j] c//? / Ji dfl = hji^ 

yrhioh l» what we wanted to prove. 


(4.6) 


It ruist bo noted that the condition postulated by this 
thoorc" can always be sutiHiied »iy a proper choice ol aigns 
for tlic loadings, it thus dotennines which oi the loadin'!;:;. 

If ®ny, should carry a negative itlt. The theta or;, however, 
docB not say anything about their nagnitudea and may be chosen 
arbitrarily to the extent that it satisiiei the lilt constrain 

•fhooraw 2* If ) and ) arc two loadings in forward 

flow satisfying the condition 2 c(u / 
then the corresponding loudin'^o 
{ % K b , 7/ ) in reverse fltnv will also satisfy 
a eirilar condition. 

I’roof: runce on<! I =f; . honce ? = und where 

the Crfly bar on toil denotes tjuantltics in revrstv flow. 
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Also 

du - da J/, „ - j IS, j A = 

'* 

Similarly . 

Henco thf* oondftlon: 

/ ( dig + ^ ! tg 

also roads 

' 2 dff / ( dfg + o/^,) ^ l( / (^ 

which provoa the t!»ef>rem. 

The tnijortanco ol this theorem is that it is sarnotimos 

i 

easier to deal wltlj n win;r, la reverse flow than in forvarrl 
flow and the above tlH;orci:i lets i;s choose the easier way. 

Theorem 3» If a loading i p, ^ ) is ropreueiitaWc by a linear 
combination of loadings (p. j I - i, 2 , .. .. , isj 

such tliat p =■ Z k P( and o< ■= Z then the 

constants are determined from the set of 
eqiiatlons Z i{ (d(k ■* ) = J(^po{^ + 

Troofs SHbstltuting the expansions for j> and « in the 

right ham! side of the above expression gives the result, 

Tlie Importance of this result is that it allows us to 
detornlne the using an integral approach which is to be 



ttfjBinsti ci ooHtjoii'tioi'i procudus"© wher© * 111 © choic© 
of the collocation polnte affect th© . The result is 
presented In tlie eanie spirit as the evaluation of Fourier 
ooifstimts in e Fourier analysis , 


Theorern 4* If th© oorroepondln'? loadin,'3:0 in revorse flow ore 



#"*0 

ip. * p. 

: °< o 

and 



. (d 

then 


and 

/■«(«/ . 

o( Cl Oc-, 

(■;. f 

Proofs 

.Since l> 

!• = P - » 

hence 

(^'k* from 


theorem 3 wo hfwo 

- J ( l>d^ * S ) dfi 

- j ( p i)!)i 4 

^ I’ p p of), * p^o()d^ 

= Z iA<.\ * ’ 

V-f 

Coraparlns; hoth sides gives the desired reault. 

As u conae<iucnce of thcorenie 3 uiiu 4 one may easily 
prove that if ( i« an, optimal loading; in fonrord 

flo^, then ( >%/ :s ) *" “I*" 

I w 


>: Colit > ‘^tij 



tt» itxpnnBlon in terns of ) hns the same coefficients 

na In tlie corrospondtn." fomwird flow, 

4.3 Final Kenarka 

The Itnportunce of the abovo results 11c In the fact, 
that tlven tm arbitrary act of lourtin'^s, It is possible to 
novo trnrards the mlninmti drn" condition with considerable 
freedom, an<l tiieroby .tiIvos the designer conaldcrublc leev/ay 
In obtalnin" near optimal wing shapes which nay be of greater 
use than the exact optlnnl in a real flow. Further, witliin 
the lliH'iirizod theory, ho hrs the option of workin”- in either 
the forward or the reverse flow depending on the particular 
probler in hand* 



ClIAPTEU 5 


MrrmonS OF IW, MTKTMIZATION-I. 

5,1 Tntrodtictlon 

In the study of opttnizerl wln^s, the problen of ohtoinins 
nlnlmura dra.'^ v/tn^a subjected to a pilven total lift seems to 
have cii!i^!{ht the v/ldest attention. This problen has generally 
been solved hy the Uer.rou^e ■"'ultlpHor method, Grahra’s 
orthO'jjoituU loadln/j; method , v.Tilch is sinnier, has fotiml 
United anpl i cation prosunnhly because the method of obtaininr; 
ortho'^oiuU l(>odljj";s os outlined by hln Is cunboroonso and 

20 

imsnitiJhle for dlT.lta] computation. The theorens of Munk 

r»/’| 0 4 

and Jones * for construe tin,"; nlnlnun drag wings have also 
found little application except for an attempt by Ginzei and 

ift 

Mnltljopp * 

ft 

Tn this chapter the l.ogrange multiplier method la 
0!itlim>8 heeiuiae of its wide use and also for the sake of 
connlete«(?aB . ’"e follow this by giving a matrix representation 

of r.raUm’a orthogonal loading method and thus make it more 
«ul table for digital computer applications. The matrix 
FD'pT'^soiititttion fiiic! oliirtod Oriih®ii^ . ' iifict his c^^worKoFs ■ sino©, 

they mention that further research is necessary to get a 

33 

li€jtt4iir olittiiiii lo-odtu^s ^ 
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Multi Iillcr Mf>t8u>r1 

The closateol Lagrari,?e multiplier (U!) laethod converts 
u conatri tned inirilnisiatlon problem to an unconstrained one, 
Tbe prototype ot this prohlon with equality constraints is 

??lninilEe 

subject to 


(5.1) 


|(£} = k . i - .w 

whorii the objective function /(’X) and the constraint 
functlone (f-J are ftinctlons of the n -vector X , The 
' ultluHer method reeulres one to form a new function 


(5.2) 


ni ~ 


(5.3) 


and to solve the following simuitanoous eouatlons 

i'P . 

(. j__ /; -M . — — = o ^ -i - 1 , ■■■■, -n 

^ 


Tdicie is the point at which /(x j takes on u relative 

mliUmtin. 


The drag mitilmlzotion problem aubjeot to a given lift 
may be stated osi 


’‘^inlnize C, 


r 2 ^ 

^ i 


(5.5) 
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subject to } €c It ~ C'l . 

r 

The corresponding La?5run«;e function Is 


(5.6) 


p ^ Cp f >t ( 2. £( L' - Ci_ ) 

i 

where the and ^ for lainimnn drag are to be found fron 

^ ^ *" ii'' ^ ''' 


(5.7) 


f'- 0 H 


C 


(5.8) 


It way also he shov/n that at the niniisum drag point 

^ / C 

(5.9) 

l.e. the Lagrange nultliilier is twice the negative value of 
the optti’i'm drag to lift ratio. 


I’ro" physical eons i derations we know is a strictly 
positive euanilty, I’urthcr, If all the constraints in the 
problori are linear, it may be shown under certain. conditions, 
to have 0 unique solution. The Lagrange nuitlplier method Is 
useful ill liticarlKOd wing theory because the drag is a cpjj dratic 
function and the constraints to wldch the wing is sJibJected 
are often stated in a linear ron!i. Thus the resulting 
oquatlons corresponding to (5.4) are linear arul may easily 
be Holvcd by any of tlio standard methods for slnultaneous 
linear equations. 

5.3 of Ortltegoiu-l Loaclings 


Grahan shov.ed that the drug in the case of orthogonal 
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loii(1tn"s v.'aH suporposahle * it was fiirther shown l>y Gmliani 
that if ono ohtoins a set of llftinr; pressure loadin,"s, o 
linear comhlnation of these, each carrying positive' lift, 
would .'Jive smaller rtran; for a ,^lven total lift then any 
Bln/?lc member oorryin?; the sane total lift, (The oorrespond- 
In/j result for non-orthoijonai loodtn'js is .'jlven by Theorem 1 
of Clini>ter 4)* It Kiay also he shown that the absolute nininun 
in dra.", may bo approaelicd ns close as desired by addinsj 
additional Itftiw; loadlnoss to the set# 

drahari’s method of ortho 'jonoU/otton 

(liven n set of loudlu'js )» ^ Graham 

proceeded to ortho-^pmnlir.e them oa follows! 

hot ^ ) bo n pressure loudin't. Two pressure 

loodln-js arc of different tyoes If 

I’i ^ ^ h ■ ( 5 . 10 ) 

vhpie e U a constant, onllod the intensity ol loodlnR. Now 

snpiiosc P-= (1?, y, , t, ) 1» 0 "ot "I no»-orthocotw.l 

louain-B end let P » ( ^ ). tn bo detomlned, 

be an ortliom.owl set derived from P , To start the process, 

Grahom chose the first member jp as 

h , h cif « ■ (5.ti) 

ff n > 

mi-' 

ITie second nenber ^ wns choncr. In tbo Inm 

k = 1’, + ‘K i 


( 5 . 12 ) 



where from the orthogonality criterion 

dfl =o 

wo have 

C - - 2 dif I ( ■'■ ) - 

Similarly the third tnenhor ^ Is obtained In the forc! 

k k " ^*4 ■" i ^ ^ 5^3 

aiifi thu orthof^onal ity oritorif^n we have 

Klsfi * =° ? 

J(k^, - ‘i'> =° -I 

or 

2 da ■+ C<^n ) ■*■ <^3 ^< 3/13 + 4/ ) --0 ^ 

(4^ 4/ ) + 4a -t % C<^23 ^ ) =0 J 

frOKi which and C^ may bo found. 


( 5 . 13 ) 


(5,14) 


(5.15) 


(5. 1C) 


( 5 . 17 ) 


Clearly; the process becomes laborious for successive 
raemhors in the set and Is not suitable for digital computers. 
Further, the orthogonal sot is rwm-tmtque since the process 
could have boon started with any raember from the non-orthogonal 
set instcod of tho first one or with any linear combinations 
of them, 

5,3,2 Orthogonal loadings ~ a matrix representation 
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Orthogonal loadings 


a matrix representation 


The dra/» cxprosalon may be written os 

Ct, ^ e'^ 


where B = [ / Is a row vector In load Inten- 


( 5 . 18 ) 


slties, and ^ ts a symmetric matrix with elements 

= i ^ dfi. 


( 5 . 10 ) 


Tile triirisf ormatlon to orthogontil loadtnjx retjulroo diC{Jionall— 
sROtlon of the matrix S . if ^ is non-singular and real, 
this li posalhle# Lot be a square matrix such that Cdi ^ ^ 
results in a dia^.onal form. lienee 



( 5 . 20 ) 

( 5 . 21 ) 

( 5 . 22 ) 


ore respectively the diagonal matrix, and the intensity vector 
in the transfomod space. 


Fror.1 matrix theory we know that B la not a unique 
tranttf ormatlon. Oneway use this non-tmiqnencss to choose 
^ to suit certain con”>utatiowil requirewonts in diagonalising 
, V,c neiition two methods of obtalnin*; , In the first 
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method, the ei.,o..vulo,.o of ^ aro olHotnod end the correspond- 
1.1?. elriMivcctors ore placed colunn liy col,.™ to constitute . 
'file dlitinnol nntrlK ^ then contains the eir.envalues In the 
dloeomil. This method Is .roll knotn end the matrix (5 In this 

ciuio 1 Iciiotvil 0 p tho ortho no rrniil nodnl rriiitrix ot ^ * ouo inoy 
Cl 1 110 rooiiiL tlmt nciy bo roclucodi to o unit nuitrix by tlie 

troruiformiition 


4' / J d:/ 


(5.23) 


'Hit* sooond riothod is to look for o matrix 6' which is 

trlan?»alur» l.o, 

^ i ‘ (5.24) 

ond the non-asoro elononts aro caicnloted frota 



whore /(• la clctomlnnnt of the {Lti) prlnolivol minor of . 
Since the «1ra.t^ is » poaltlvo definite? fnnctlon, none of the 

are zero. The second method Inplies that the first column 
vectoi of c;' is nsaiined en incident with the first louUtn,", 
the aoconti column vector is chosen to lie in the plane detcr- 
nlned hy the first end second load liie;, the third column vector 
is confined within the suhsoace determined by the first, 
second and third loadlnfsa, and so on. one rooliaos that this 
is what the orl,"inal Orahoiti’s nctliod was tryinf^to do. ’ , . 



A twii coiis fci'ii £ nirrizut io n n3rol)Jci":i 


5;3*3 

It may now bo, expected' that the drafi; minimization 
C0.,lcn 1 II 1 1 0 Ji wot! I cl be sij’ipj Lfiocl* Consider ^ for cxainplcj tfie 
pi 1*0 1)1 ci'*i of o b t a I li i„ I \'\ i 1 1 1 m i ,d ra ^ s \ ib j oc t to a '**; i ve n 1 i f t b uci 
pltchin/; tiiomcoit constraint* The problem is 


Mini mile 

(d ~ 

I 4 <iu 

i 

siihjcot to 




Z €i 

k 1 


i 



z 

t 

■mi =* L,»t j 

where du , 

t, 

are the timo;, lil't and pttchimr nmnent 



• ih 


coef n c I ciitH i*OHpocti vc.ly dtie to- nni t . intensi t.y of tlie 


mei:i1»cr in l.lu‘ ortlio'tntia] set. nnd nrc Urn nrescr H km! 


lift and pitchin'’: monent coclTtcients 
correspnrid I ri'!; !inc<inHtra i ikmI proT>lei»i is 


Minimize 


leNa H- a, fj + t5.:iT) 


The nccossory caadi ti 'iiis (5.4) show that 

2 du +• t =^j i= i/, ■■■■ , 'n. 

1 . 

From (5. 2B) and (5.20) one may ohtuin 
2Cj, -f- " 


( 5 . 28 } 

(5.20) 


O 
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or 




( \ t m ^ ^ j / 


SubstltuttW! for In (s.ag) no obtain 




c ^ 


(5,30) 


(5.31) 


and sulmtltutlns for fe. in (6.29) we obtain 


c; 


d / ^ i- 4- 

7'% f J °'" 


j 1 I "i - / z 

/• ^ii j ij ^ •'jj 




h' \ ^ 

1) 


(5.32) 


and solvin'^, for A, glvoe 


^ C'ft) 

a 


T 


i’lf 


^ J i 5 

i <ij 


z 


1 


(5.33) 


Ecjuatione (5*31) to (6,33) give the winimun value of the drag 
anr! the required load Intenoltioe , setting =o givee 
the reuni tu obtained In reforenoo 32 for the cobc when only 
the lift oonatrnlnt la present. 


Ab a finol remnrlc, the Lagrange milttplier method xvhen 
applied to non-orthogonta loodin'»H may give riue to an ill- 
coridltloiad matrix and higher precision may be required to 
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to ohtniri roasouablo rosuits. 'rhe ortho.";oiio;i. Inodim' Rjetiiod 
alloviatoB tliis problem him! should be more suitable for 
nuraerical claculatlons . An additional advantage here is that 
one nay el\irnys choose a proper sign for ^ so that ^ 
posltlt'o, lioiice the intensities according to Graham must 
all b(' non-ncgntlvc for the riinimun drag problem subjoct to 
a ffivoii lift. This is u useful feature in nathenntical 
programming formulations ns will he shown in the next chapter. 
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CIIAPTKR 6 

FItAGTICAl. }!I TIK)DS OP DI5AG MIKIMIZATION-II 


G*l' Introduction 

Fo continue in this chtipter* the study of dm , *5 minini— 
zatlon teciuiiqiies, and explore the poosihUity of using 
r-atliw^'otlcal progmnnlng (MP) in order to Introduco inequality 
oonstrntnto in odditlon to the equollty constraints comwonly 
used,* Mnthrnntloia progmnrilng has found wide applloations 
la t!ic optinlzatton of strnctiircsi hut very little in the nero- 
ilynj riles of wings and bodies, ny using MP one nay Introduce 
such constraints In the design problem ns, for example, a 
snitcihle twist distribution near the wing tip for favourable 
tip stall I etc* 

* 

The general tiP problem Is 


Mliiiralzo 
subject to 




and 



4 

i =12, w 


% it) 

^ hi 

i 2 UO, - • 7 F- 


% iv 

- k 

J . lH), - ■ . ”> 

( 6 . 1 ) 


^ 0 ■ ■ 

J " ■ 



reviewer has brought to the author* s notice, after stibmisslor 
of this dissertation, that MP has been used to solve some thick- 
ness {.rohlons In wim-horty nnrodjmnnlcs In n 

.7. ABtronnntleol Soi., nn.283-296, lOCB. Ko hove used MP to 

solve lifting wing problems. 
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A 1 , 0.1 rontun- of this formnlatlon Is that the tmknonm 

vorlohloa Xj «ro ronolrod to be non-ncf-atlvc. Thte allows 
Oise to obtoiis i.l f^orl thr'B tlust woric more efficiently by res-* 
trl (5 1 1 iif.r, the search to one orthant* Iloweverf any variable 


with unrestricted stfsn may be put in this forra by definln;? 
tivo now iioin-iio,«;otivo variables x' , x", so that 




X.: 


J ’ J 


■Xj 

z- o . 


( 6 . 2 ) 


A prohlof” 1)1 wbtob i(l) and all the are iln(*ur is 

called a linear pro":r£itnnin" (LP) problon. A problcn in which 
at least one of these functions is non-llnoar or has spoclol 
rostrictinnu plnecrt on it (o./;, sono or all Xj nicy bo required 
to he Inteners) is called a non*lincnr prosrarinjing (NhP) 
problew. A special case when -f^l) Is a quadratic function and 
nil the ore linear is called o qimdrs.tic program lng(Ql*) 

prohlen. The draw; mlntmlsRatlon problem subjeot to linear 
oonstrnlnts falls under the QP formulation. However, approxi- 
mate LP fornulntlons are also possible. 


6#2 Linear hronrawming: MethM 

The problem under consideration is restricted to the 
conical wiir'j problem and for wln/^s in subsonic flow. This 
problem allows control over the Induced drag which Is a ntiiiiruaa 
If the opanwiHC dlotrlbiition of lift is oiliptic. The fornu- 
Intion, therefore, reenlren u act of loa.lin't.s to adjust tin tr 



87 


intenslti <?8 such thnt It satisfies the constraints In the 
proMcn, onrl follows the Ideal oUiptlc epanwlse distribution 
as closely as possible* 


As an oxonplc, the prohlon of nlnlinum drag subject to 
a given total lift Is 


MininilsiQ 
such that 


(r>.3) 


^ (Ji ) *“ ’2 ^ ^ ^ 

£(t; + Z'S; CjLD s S J 


(0.4) 




c. 


(0.6) 


^ , o 


((i.6) 


Where (0,4) signify, in a discrotised fortn, the constraint 

j£(l) - CjCD j 6 ?. (6.T) 

and are n-spanrlio oolloootlon points. E('^ ) is the 
©Hlptle load tliatrihutlon and Cj('7) i® the spanwlse load 
distribution for the J**" loading of unit intensity. For the 
best fit obviously S tnnst bo ratnlnlsed. Also nottoc that 
. This will be so If the chosen loadings fona an 
orthogonal sot with each loading carrying positive lift, ror 
other oases one nay either orthogonallzc the sot or use ( 0 . 2 ) . 

The solution to the prohler. will give the intensities 
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and the laaxinura deviation S between the Ideal elliptic dis— 
trlhiitlon and the actual dintributlon* The smal lex* the value 
of 6 I Bpaller will be the drag which nay be calculated once 
the Cj ajre known* In the next chapter* will be shown that 
this Ricthod Is good. 

Additional equality and inequality constraints are cosy 
to InoorjKirate into the problem. These constraints nay be 
usodp for example I to restrict the solution to regions where 
the theoretlool oeaumptions are valid, (me night also require t 
for instnitoei the angle of attack not exceed a oertain value* 
^Cti^ f i • e • 

I ^ I - • ( 6 , 8 ) 

Till* may be treated as two Inequality constraints 

i 

One might further introduce the constraint that should 

progreesivcly dccrenec os one moves from the wing root to the 
tip, which In a dlsorcttacd form would bo 

r«.- “i ^ V ^ ° , (o.Ki) 

I 

where >l^jl . Such a twist distribution sliould al levin te 

problems of tip stall. 




80 


Viiirlous other constraints orislnf^ due to inanufacturliig 
roqiilronicnts may he constructed and added to the constraints 
list in the prohlctn* Thus the method linnodiatcly opens up a 
wldor possihUity of obtaining roollstio wing shapes. 


Several methods exist for solving Li problems. The 
output of whioli liKUeatos whether the solution is feasible or 
not I if yes# whether it is bounded or not* Xf the solution 
io feoslblo end bounded i it gives the values of and S* 



The firm being a ijiiadraUo function of the Intensities 
, and within the frantwoilc of the Unearistd theory, many 
praetienl eonstralnts ean be expressed linearly in terms of 
£■ , In this farm ths problem is 

i 

UiniMisa Q> s ^ ^ '^y <d.ii) 

» j 

subject t« 

r ^ (0.12) 

and other linear equality and inequality eonstralnts some of 
wliieli may be (6.0) and (0.10). 

If one uses orthogonal loadings with each member oarrylng 
a positive lift, the problem may be written ns 


Itlnlmlxe 


( 6 . 13 ) 
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subject to 


= £ • : 



i 


(6.14) 


4 ^.. , t 


(6,15) 

whore (^.15) 

mpmBBBB 

the additional linear constraints 

on 

the problem, 

it aiqr. 

The non-negativity condition on 

by 


virtue of ortho .'.^omaity provided each /^ > c? ts nutomotically 
sail of led* 

Sovorol methods for solvlnn; Qi» prohlons exist in which 
the ohjjcctlvo fimotion is the sum of a linear and a quuUrutlo 
function of the unlcnown variables. I’rorainent nnon," these are 
the mothodH of Beale, and Wolfe which may be used. 

6.4 Hemcrlcs 

A 

The Lan/ironr^o multiplier method has been widely used in 
drag minimization problems. In comijarison, the application of 
pro'^rnmmln" methods Is of recent ori5:in. It was shown, in 
this context, how hP and QP formnlnttons may be usefully 
employed In dra'? minimization problems under Inequality 
constraints and how orthogonal loadin'^s may be useful In 
reducing the oomriloxlty of the j*roblera. Use of orthoRonal 
loadin/i^s also gives substantial benefits In computer core 
rcciutrcwonts* 

IThoreas the QP formulations are straightforward, the LP 
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foiTiulotloiifsi require one to choose a suitahlo objective 
funotlori so that the resulting problem Implies u drag mini- 
miantion problem. The elliptic spanwise distribution of lift 
for mlnimuKi induced drag hos been used in this chapter. 
Ilowevori on similar lines other drag mininizotlon theorems 
may also be used* o.g, the combined flow field results of 
r.unic nnd Jones,' This olso provides a convenient vehicle for 
(lotervilnlng how sensitive nlnimun dreg shopes ore to smell 
excursions from the theoretlonl mintimim. 



92 


CHAPTER 7 

COMPARISON OP 0RAO MINIMIZATION METHODS 
AND OPTIMAL CONICAL WINGS 


7*1 Introduotlon 

A lew examples are ehoaen to llluatrate the drag ralni- 
wlsaiion nethoda eutllnei! in chapters 5 and 6. Also studied 
In this eh«M;>ter are the wing shapes carrying a given total 
lift usini't vArloiis oonhinations of the first four haslc twist 
dlitrlbutleiie in the oonlosl series of chapter 2 to provide 
daaign infortoation* All examples ore for <a»o€* and tan^-Zo, 


Comp aratt ve Study throogh 


rihcaiaple 1. Miniaiixe the drag subject to a given ^o l » 
using the first three haslo conical twist distributions, i.e. 


Minimise 




-ii 


V 



subject to 

1 1 *. /.• = ' ' 

Where Is related to the £;jri through Equation (2*25) 

with 'The quantities dgj , and are oelculated from 

Kquationa (2.40)-(2.43) . 


After substituting all the numbers we have, 
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MinitnJiseo 


- G-mHs €f - 0-n 0-171) _ o-o72/s6 £ €, 

“ 0-OZL,iio ^ - o ozasiz Cif - O-oOA-JSf £4 

* 


tubjeot tf» 


0 ■!• Q-(*Z8l^b £^ t 0-2765iS £fc = 0-1 


W« have alraady noted in chapter 6 that mathematical programm- 
ing methode require an additional constraint that all the 
unknown variahlea he non-negative. If this condition cannot 
be guaranteed apriorii an equivalent problem nay be described 
using Kqwatlon (6.2) to define a set of new aon>negative 
variables. An alternative method is to convert the given 
set of loadings to on orthogonal set using any of the methods 
given in chapter S. Here we ohoose the di agonal i*at ion method 
using the eigenvaldUe technique. The above problem is now 


rewritten as 
Minimi see Cp 


Iklih 0-3SS9c>x<0 


subject to 




0- « 6737 xft' €i - - 


The matrix 6 of Kquatlon (5.20) that accomplishes this is 

” 0.807299 -4l,4836if 0,117936 

G ■ 0.423886 0,593294 -0.684373 

0.260988 0,643521 0.719561 
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.. * . ' 

Tli« aolutlon to this problem is 

/ ' t 

« 0,0897, 0,4331, 3,1712 

and 

Cp m 0,001294* 

This wa« aolved ualng Qp and the UI methods and was used to 
oheok out our Qi’ computer program. 

However* to eolv# the above problem ualng the LP method 
requires a linear objeotive function to replace the quadratic 
drag function. For this we follow the method outlined in 
chapter 6, fbus we write 

Minimise € 

subject to 

-c'c^j) (i - oijj (<! ■ 

for j= a. •■ - ^ 

and ^ j , 

£ €« , ^4 - 3 ^ ^ • 

Her© € is a positive quantity to be minimised, ^(-V) is the 
elliptic spanwise loading - a condition for minimiw dreg, 

C^l Cl) are the spanwise loadings due to the orthogonal 
loading of unit intensity and are the collocation 
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pout. «l.oro tho o.Mltlo. (6.T) U to bo satisfied. 

several value, of n .ere cho.en. A sample of the 
results is given below: 


1) * 3 : * 0.0 » 0,3 

* 0,0890 1 0,5132, 

Cp » 0.001304 

2 ) ■H m 6 ; '/?, « 0 . 0 , 0.1 

. I , 

a 0.0874, i^m 0.4800, 
Cp m 0,001302 


0.8 

t 

S « 0.000246 


0.3, 0.4, 0.8, 0.9 

/ 

4* 4.5451, S s 0.004223 


3) ;i s» 10 J m 0.0 to 0.9 In steps of 0.1 

■ 0.0877, m 0.4697, C^xx 4,4582, 0,004534 

Cp m 0.001300 . 

Coraparlm?: the results with those of the QP (or LM) method we 
find that with Increasing /t , the Cp approaches the correct 
mtnlftuiwi. It may, therefore, be concluded that the LP method 
Is a useful technique, 

7,3 Mtntmuin Drag Wln^gs with Inequality Constraints 

Minlwuni drag shapes under inequality constraints are of 
groat practical significance, A few examples are given to 
illustrate the QP method. 

Example 1, Oslng the basic twist distributions 1 and 3 
wlnlmlase tho drag subject to a given o f , a root incidence 
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between 0.03 and o,04 radinn 

v',u« raaian, and a continuous washout 

in incidence from the wing root to the tip. 

After orthogonallaing the basic loads 1 and 3 end 
dlsoretlrtng the washout constraint, the problem *oy be 
fltatetfs 

|3L, ^ 3i. 

Minimice Cp '-s o ili27 i 10 dj 

subject to 

€* -h O' » / U I A /o"' » 01 




44#! 


^ O 


■i = C, 2 , .... , H 


£ii{{f 4 O' 0 4 


where 




0(l ^ Oi (rl. ) 


a o- a 


.fj 

114 f 

1,0 

% 

The results are 


C Of. . 0, _ 

~ 0 • 3^^ 

tt f-O 


ih -> » 


* O.ilOiO, 0.15443, Cp a 0,(»01418 

In ooiiipurl.nn o wlni? designed only tor would have 

fj, = 1) ,001321. Thu. the additional eonstralnts In the exanplo 
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have increaafd the drag coefficient by about H percent, 
tn a real vleoous medium thia win? Is likely to be better 
than the one designed with only the lift constraint which 
gives a wavy shape and a washout only near the tip. The 
wing shn;ie and Incidence distribution ore shown in Pig, 7.1 
and compared against the case with only the lift constraint, 

This example shows that it Is possible with only two 
hcslo conical shapes to design a wing which satisfies a 
dsslgner's notion of the desirable properties of a good wing 
in a viscous medium even by using an Invisold analysis* if 
one introduces inequality constraints. If ons hod solved 
the classical prohlsm of drag ratnimiaation with only the lift 
ocnstralni» this prcetioally useful wing would not have been 
indicated# 

Example 2, »sing the basic twist distributions 3, 3 and 4 

wlnlmiae the drag fer n given CwO.lt wing root incidenoe 

radian, 4 0 from wing root to tip and « o 

for 0,8 ®hd Z if' *»0,04 at '») « 1,0. 

Orthogonalialiig the basic loads and dlscrettxing the 
prnbelm wo have the follovfings 

Cp * <?• > 0-8lSU6ttn % ^ 

-6 ,2 


Mini mix© 






NS OF EXAMPLES OF 1 . 82 . 
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siibjeot to 


where 


Z- 

- /i 


•" o/j- ■£: o 

■Vo 

= 0-0 3 

ni/l 

- 0 - 0 9 

i 



^4, S ^ 


(■ - o, i ■■ ■ , io 


5^ io 


<y. . o/l-rj- ; _ 


and the ^t' as in Example 1 above. 


The resulta are 

0.17389, 0.40491, 17.05696 

<5j S 0.001317 

But with only the lift constraint, Cp = 0.001277. 

As for the previous example, the results are shown In 
Fig. 7.1. In this example the wing is restricted to have 

negative spanwise surface slope, and for favourable tip stall 
characteristics, a washout is prescribed in the outboard 
portion. The last constraint Zff ^ -0.04 in the example 

means that a minimum clearance from a datum for the tip is 
specified . 

This example again shows that our basic shapes are 
capable of configuring wing shapes to practical requirements 
with few members. Such problems for lifting wings do not 



too 

appear to hove been solved earlier in the available literature. 
Never thel ess j these are the problems that are faoed by wing 
designers, and Mf methods have thi potentiality of providing 
the necessary solutions. 

Tnctdentolly in the above examples, the LP method could 
also have been used as indicated in chapter 6 instead of the 
more accurate QP method used here. 

7.4 Mlnimtim Drag Wings for Given Lift 

Wo close this chapter by studying the classical problem 
of minimum drag shapes obtainable by superimposing two or more 
of the first four basic distributions of chapter 2 in all 
possible combinations for a given total lift. Many of these 
shapes appear to be of practical use. The results are 
summarized in Tables (7,i) end (7.2). The wing sliopes are 
drawn in Pig. 7 *2 . 

Very systematic trends are observed as combinations of 
two, three and four of the basic twist distributions are used. 
For on even number of combinations it is noticed that the wing 
ordinates He above the reference '^j -axis and for odd number 
of combinations they lie below.* Combinations of higher order 


* Combinations of greater than four basic shapes were also 
calculated but are not shown in Tobies (7.1) and (7.2). 
However, the comments made in this section reflect these 
resul ts also . 
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TABLE 7.1 s AERODYNAMIC AND GEOMITRIC EROPEimES 
OF OPTIMAL CONICAL WINGS (a « 0.5) 


Wlnga* 

12 

13 

14 

23 

24 

34 

V 







0.0000 

0.6381 

0,6389 

0*6395 

0.6395 

0.6396 

0,6396 

0.1000 

0,6359 

0.6388 

0*6404 

0.6413 

0.6419 

0.6422 

o.aooo 

0 # ftS44 

0.6405 

0*6442 

0.6462 

0,6483 

0.6503 

0i3000 

0*6492 

0*6516 

0*6558 

0.6535 

0,6583 

0*6628 

0.4000 

0,7032 

0.6870 

0*6850 

0.6706 

0*6753 

0,6796 

0*5000 

0*0224 

0.7718 

0.7508 

0.T216 

0.7152 

0,7084 

0*6000 

1 .0368 

0.9400 

0,8874 

0.8544 

0,8189 

0.7835 

0*7000 

1.3148 

1.2234 

1,1450 

1,1335 

1,0622 

0.9916 

0,8000 

1 .6276 

1.6098 

1.5613 

1.5929 

1,5300 

1,4683 

0,9000 

1.7485 

1.8991 

1.9969 

2 ,0485 

2,1203 

2,1924 

O.OJHM) 

1,5186 

1.7417 

1.9312 

1.9625 

2,1349 

2,3065 

0,9060 

0.5760 

0.6956 

0 .8135 

0,8139 

0.9302 

1,0464 

1,0000 

0.0000 

0.0000 

0 .0000 

0.0000 

0.0000 

0 .0000 


ifc stand* for wln« daaiened to liava a 

, ■ . a(, aifc 

twlat dlstrtDutton^ * % "? t * * , 

e.g. Wing 12 has - v 
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TABLE 7*1 (Conta.) 


i~cn) jc^ 


Wlttgii 

V 

123 

124 

134 

234 

1234 

Flat 

Plate 

0*0000 

0,6829 

0.6785 

0.6749 

0.6710 

0.6460 

0.6368 

0,1000 

0,6918 

0.6864 

0.6816 

0.6761 

0.6458 

0.6398 

0.2000 

0.7108 

0.7051 

0.7001 

0.6939 

0.6570 

0.6497 

0.3000 

0.7190 

0.7195 

0.7215 

0.7239 

0.6995 

0.6674 

0.4000 

0.7016 

0.7144 

0.7305 

0.7523 

0 .7692 

0.6946 

0,5000 

0.6735 

0,6924 

0.7162 

0.7526 

0.8189 

0.7351 

0.6000 

0.7046 

0.7024 

0.7050 

0.7227 

0.7957 

0.7958 

0.7000 

0.9237 

0.8693 

0.8156 

0.7688 

0.7655 

0.8914 

0.8000 

1.4642 

1,3797 

1.2888 

1 .1864 

1.0774 

1 .0610 

0.9000 

2.2030 

2.2343 

2.2660 

2.2546 

2.2205 

1.4605 

0.9500 

2.2663 

2.4164 

2.5616 

2 . 6880 

2.8085 

2.0388 

0.9950 

0.9938 

1.1135 

1.2311 

1.3461 

1.4761 

6.3741 

1 .0000 

0,0000 

0.0000 

0,0000 

0.0000 

0.0000 
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TAWJ-; 7.1 (Contfl.) 

* 


Z fanG/ % Ci_ 


wings 

1 

12 

13 

14 

23 

24 

34 

0.0000 

0 .0004 

0.0002 

0.0001 

0.0000 

0.0000 

0 .0000 

oaooo 

0.0088 

0.0044 

0.0029 

0.0000 

0,0000 

0.0000 

0.2000 

0 .0346 

0.0175 

0,0115 

0,0007 

0.0004 

0.0000 

0 . 3000 

0,0768 

0.0393 

0.0260 

0.0035 

0.0019 

0.0003 

0,4000 

0.1293 

0.0692 

0.0461 

0 .0105 

0.0059 

0.0014 

0 . 5000 

0.1911 

0 . 1062 

0.0716 

0.0231 

0.0139 

0.0049 

0,6000 

0,2560 

0.1478 

0.1016 

0.0419 

0.0271 

0.0126 

0.7000 

0.3175 

0.1895 

0.1334 

0.0644 

0.044" 

0.0255 

0,8000 

0.3682 

0.2244 

0.1614 

0.0838 

0.0619 

0.0406 

0.9000 

0.3993 

0.2412 

0.1739 

0.0868 

0.0657 

0.0452 

0,9500 

0.4044 

0.2381 

0.1682 

0.0758 

0 .0548 

0.0346 

0.9950 

0.4016 

0.2261 

0.1521 

0.0547 

0.0323 

0.0107 

1.0000 

0.3984 

0.2241 

C .1496 

0.0516 

0,0289 

0,0060 
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TABLK 7 a ( Contd .) 


I Bn 9 j^<^L 


Wlilg » 

123 

124 

134 

234 

1234 

Plat 

Plate 

0.0000 

- 0.0006 

- 0,0004 

- 0.0002 

- 0.0000 

- 0.0005 

0.0000 

0 , 1000 

- 0.0120 

- 0.0087 

- 0.0050 

- 0.0001 

- 0.0119 

0.0000 

0.2000 

- 0.0462 

- 0.0338 

- 0 .0200 

- 0.0015 

0.0443 

o.oooo 

0.3000 

- 0.0976 

- 0.0723 

- 0 .0446 

- 0,0068 

0.0883 

0.0000 

0*4000 

- 0.1681 

- 0.1197 

- 0.0770 

- 0.0185 

0 . 1335 

0.0000 

0*5000 

- 0.2108 

- 0.1696 

- 0,1139 

- 0 .0365 

0.1721 

0.0000 

0.6000 

- 0.2753 

- 0.2157 

- 0,1497 

- 0.0567 

0.2035 

0 .0000 

0.7000 

- 0.3216 

- 0.2537 

- 0.1786 

- 0.0722 

0.2348 

m 

0 .0000 

0.8000 

- 0.3632 

- 0.2852 

- 0.1992 

- 0.0780 

0,2753 

0.0000 

0.0000 

- 0.4162 

- 0.3244 

- 0.2238 

- 0.0838 

0.3205 

0.0000 

0.0500 

- 0.4563 

- 0.3573 

- 0.2490 

- 0 .0988 

0.3323 

0 .0000 

0.0960 

- 0 , 6062 

- 0.4023 

- 0.2885 

- 0.1306 

0.3232 

0 .0000 

1 .0000 

- 0.5129 

- 0.4086 

- 0.2943 

- 0.1429 

0.3204 

0 .0000 
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TAOLK 7.1 (Contd.) 




U ton 





wing. 

1 

12 

13 

14 

23 

24 

34 

0.0000 

0,1568 

0.1698 

0.1764 

0.1824 

0.1855 

0,1885 

0.1000 

0.1655 

0.1742 

0.1793 

0.1826 

0.1856 

0.1885 

0.2000 

0,1899 

0.1873 

0.1879 

0,1845 

0.1866 

0.1886 

0.3000 

0,2247 

0.2085 

0.2023 

0.1924 

0.1911 

0.1698 

0,4000 

0,2612 

0.23S5 

0.2219 

0.2102 

0.2027 

0.1951 

0,5000 

0.2873 

0.2626 

0.2447 

0.2385 

0.2241 

0,2099 

0.6000 

0.2871 

0.2776 

0.2639 

0.2685 

0.2529 

0.2377 

0,7000 

0.2416 

0.2586 

0.2617 

0.2756 

0.2719 

It 

0.2685 

0.8000 

0.1280 

0.1697 

0.1978 

0.2110 

0.2323 

0.2536 

0.0000 

-0.0798 

-0.0440 

-0.0089 

-0 .0087 

0.0260 

0 .0604 

0.0500 

-0.2282 

-0.2215 

-0.2089 

-0,2149 

-0.1996 

-0-1845 

0.9550 

-0.3915 

-0.4349 

-0.4730 

-0.4779 

-0.5133 

-0.5484 

1 ,0000 

-0.4210 

-0.4622 

-0.5084 

-0.5124 

-0.5562 

-0 , 5996 
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TABLE Ta (Contd.) 




oC tan 

S /<i. 




Winga 

1 

123 

124 

134 

234 

1234 

Flat 

Plate 

0.0000 

0.2485 

0.3409 

0.2330 

0.2225 

0.1802 

0,1927 

0,1000 

0,2368 

0.2324 

0.2282 

0.2222 

0.1911 

0 .1927 

0,2000 

0.20T1 

0.2097 

0.2134 

0.2183 

0.2154 

0,1927 

0,3000 

0.1T39 

0.1817 

0,1910 

0.2046 

0.2303 

0.1927 

0 . 4000 

0.1565 

0.1617 

0.1665 

0,1805 

0.2170 

0.1927 

0.5000 

O.ITOT 

0.1653 

0.1603 

0.1585 

0.1796 

0,1927 

0.6000 

0,2184 

0.2018 

0.1839 

0.164C 

0.1540 

0,1927 

O.TOOO 

0.2734 

0.2590 

0.2419 

0.2191 

0.1883 

0 . 1927 

0,8000 

0.2643 

0.2743 

0.2320 

0.2838 

o!2730 

0.1927 

0,9000 

0.0542 

0.0887 

0.1233 

0.1570 

0,1910 

0.1927 

0.9600 

-0,1925 

-0.1722 

-0.1509 

-0.1274 

-0 .0986 

0.1927 

0.9950 

-0.5356 

-0 . 5680 

-0.5993 

-0.6278 

-0.6579 

0.1927 

1,0000 

-0.5823 

-0.6239 

—0 .6646 

-0.7105 

-0.7443 

0.1027 
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TABLE 7. I (Contd*) 




tz 13 14 33 24 34 


Wlii«?a 

V 

0,0000 
oaooo 
0.2000 
0.3000 
0.4000 
0.5000 
0.6000 
O.TOOO 
0.8000 
0.0000 
0.9500 
0.9960 
1 .9000 


0.6361 
0,6408 
0.6468 
0.6500 
0.6447 
0.6179 
0.5593 
0,4596 
0.3157 
0.1402 
0 .0555 
0,0019 
0.0000 


0,6389 

0.6392 

0.6383 

0.6369 

0.6276 

0.6043 

0.5567 

0.4721 

0.3391 

0,1594 

0.0651 

0.0023 

0.9000 


0.6395 
0.6387 
0.6359 
0.6297 
0.6173 
0.5938 
0.5511 
0.4767 
0.3547 
0.1754 
0.0739 
0,0027 
0 .0000 


0.6395 
0.6373 
0.6317 
0 . 6230 
0.6107 
0.5909 
0.5542 
0.4846 
0.3624 
0.1784 
0.0746 
0,0027 
0.9000 


0,6396 

0.6371 

0.6302 

0.0191 

0 . 6038 

0.5821 

0.5473 

0.4854 

0.3740 

0.1928 

0.0830 

0.0031 

0.0000 


0.6397 
0.6365 
0.6293 
0.6160 
0.5972 
0. 57 25 
0 . 5407 
G.4S73 
0.3854 
0.2067 
0.0924 
0.0032 
0.0000 
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TABLE 1A (Contd,) 


‘eCv) 


123 


124 


134 


234 


1234 


Flat 

Plate 


Wings 

'i 

0*0000 
0 . 1000 
0,2000 
0*3000 
0.4000 
0,5000 
0.6000 
0.7000 
0.8000 
0.9000 
0,9500 
0.9550 
1.0000 


0.6829 

0.0737 

0,6462 

0.6113 

0.5774 

0.5498 

0,5240 

0.4803 

0.3843 

0.3039 

0.886 

0.0034 

0.0000 


0.6785 
0.6697 
0,6461 
0,6129 
0.5775 
0 . 5457 
0.5167 
0.4760 
0.3901 
0.2162 
0 .0965 
0.0038 
0.0000 


0,6749 
0.6675 
0 . 6470 
0.6157 
0.5785 
0.5417 
0.5085 
0,4701 
0,3944 
0,2278 
0.1044 
0 .0042 
0 .0000 


0.6710 

0.6657 

0.6496 

0.6214 

0 . 5 S 26 

0.5302 

0.4993 

0.4610 

0.3953 

0.2381 

0,1118 

0.0046 

0.0000 


0.6461 
0 . 6469 
0.6455 
0.6321 
0 . 599 5 
0.5502 
0.4973 
0.4514 
0.3926 
0.247 5 
0.1196 
0.0050 
0 .0000 


0.6368 
0.6334 
0 . 62 38 
0.6073 
0.5835 
0.5513 
0.5093 
0.4546 
0.3820 
0.277 5 
0.1088 
0.0636 
0,0000 



Flat 00124 80.76 1.5560 ,01027 -0.01027 

Plate 
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twlBt dl.trlh..tlo«. in either case, tends to flatten the sine, 
surface earlno.s, apart fro,, the leadln* edge droop, although 
dtsocrnlhle. Is saall. The n»her of waves Increases with 
increasing nunhor of hasic twist distributions. All the 
Shape, have leading edge droop which should ^ -2:; 

flow, and the outboard portions have a washout in no 

which .held alleviate tip stall. 

Table (T.2) *“"* ^^It^wlth 

.e.olents culte close to the plane delta with 
obtain drag coefUolents quit 

in 7.2 show that tneso 

suction and the optimal s op 

haoes are practically useful In the sense that they 

shapes are pre« viscous flow. 

V + Inna of a Rood design in 
vloiat® our notl i 

, Tsble <7.1) fat uae of higher order basic 

It is clear twtata in a design 

l"”e.sure pesh outboard. Increases the peaK tntenstty 

ling on the forward facing drooped portion of the 

and hy acting * r nn that alleviates drag. 

a.inPB a thrust force that aiie 

loading edge produce in«»t bv 

ieapj. 1 % ^ _ gwctlon force lost oy 

in fact Table ( .^) 

pneecrlhtng an .„.d 

hly mode up the plane delta wing with suction 

arc attaimdJl®' 

interesting feature of these designs Is that 
Another Interestl the wing Incidence 

ivon of twist distributions xn« 

for an even number 
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and thcf wlnA root pressure coeXflcient increases as higher 
order hnslo shapes are used. For odd number of corabinatlons 
the opposite trend Is true, Table (7.1). Further, unlike 
the works of Tsien and Holla, where marked variations In the 
optimal shapes may occur with changes In the design Mach 
niimher (or the slenderness parameter) the present designs 
show rather mild changes (data not presented) . This may be 
taken to indicate that the off-design performance of our 
wings will he hotter* in the sense, that they do not show 
abrupt changes in behaviour, and hence are likely to be more 

useful In practice. 
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CHAPTEB 8 


discussions and extensions 

1,1 ■lYpr*"' Pr«»«nt 

8«iror«l problem oonneoted with tb« doolsn ■>* 

, .1 oai»b«r hov. bom .tudled. Design nethods 
«ing using oontoal oonber ,„„„n»iity 

HSV. been studied u.. 

oomtrulnts tor obtslnlng ■tnl’"” 

Illustrated, tor a conically conbcred delta .lag 

.d«. in a supersonic strean. and carrying 
.nbsonlo leading edges In 

, , -t.t distribution, escb term ot ubich an 
a polynoulal tulst d additional 

in the eonloal coordinate. An 

even orde leading edge is sero at 

^ A * 4 «« In that the lo«« 

„.trl«tlon ^ ,„etlc.ny usetul condition 

♦hil design condition. 

the des g ^ aeparatlon.and the 

gin©# gttoti ^ real 

.-e-iv to be closer to the real 

theoretical results are ontalaed 

chls way useful wing shapes have been 

aitufttlon. ftioae to the plana 

A characteristics come close xo w k 
whose aerodynamic chara 

acitn with leading edge suction. 

* A and 6 are, however, not res- 
mtsb results ot Chapters 4 and 

fhe resuiv» -«nlled to more general 

irvwiiv but may 

trloted to delta w ngs . gubsonlo and superaonic 

.ad twist distributions In subsonic 

rosnlt. ot tbo leading edge vortex pbenomenon In 

aroni.lng and with further refinement. 

Chapter 3 appear promising 
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mny torn o roUnMo nothod of predlctln.. proBsureo, «t leost 
for the 0...0 of o plone dclto wl.« ot on Inoidonco. It olBO 

InciiKloa conprosailJiilty crxocts. 

In Chapter 7, the pro^ramnin,'; methods arc illustrated 
hy nur^rteal examples lor nLni..u.n dra. win. designs subjected 
til realietlc oimstralnts. They are thus shovn to he 
«leMl*'n tools. 


8.2 *..itftpe hxtensteM 


,t tmuld he of intereot to obtoln the lo« sneod ohoroo- 

*1 itAnioal wln-’B described In Chapter 2 and to 

torts tics oT the conical 

4 - Mild of nodlflcatlons they ml'dit need, if any. 

800 frtmt hlnn oi noui* 

«t ratios, the leading ed.e vortex phenomenon 1. to 
1„„ nepoct rotloB. ^ ^ ,,,trlhnt.on nUl 

o...ectcd. ft fB not yot 

nffeet the — ^ f„ .fn oron 

With . thlcuucss d S the phe«<n.enon 

48 definitely needed, not oniy 

lU CO! I'i riaoit v;lth 

„ Olc- rly. hut nlBO to -.onerntc doW 
I'jorc cie* r 1} • 

theory. 

VA ed. Li* CUci tor d hr'vc 

.,r,.-,ru™,ln tochninnoB PfO.o.cd 

.. to nrohlc B ,vi th no.-lincor const,. xntB. 

n notnroj <r,' ^ 

iicc'u J * 

|..,rthor the ohjee IVO ‘ e«.' .1<-, It nny ho 

.s+iif'f critorii . i o* ■ 

4 .h Kill* iiltfiljC OtiHCI 

oxproBB.o., ^ eonhlnotlon of .Irnh 

the nonMfnotnrlnf • ^ cortnln rofiht 

„„„ ..cnnfaotnrlnt cost «IU. oaoh 
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fnotor whloli OTUli) Imllootc tholr reli tlvo Inportnnoc in the 
D'.tlnlnotlon aoolnions. 

Sevcrol toohnlqnes nre knonn for solvlnf; the rjenernl 
ron-Jliiwr iro’rnnnln", vrohlc" whore the ohjectlvo fenotion 
e„.1 the ooiwtrolnts rre nllowocl to ho non-linear, locccpt 
nmlcr rentriotlvc olrcnnntonces, none oi thene r..norontDO a 
-lolna optl. « nointion. They ;;cnernliy -,ive .n loom o-.tlro.. 
whleh nay or nuy hot he the mohnl optl.e.r.. In any e..oe. it 

-noro than likely that the eolntion wool,, he ooeTol or a 

Tnr the llrst Iterutton wlli 

iltfloront «t the snltrtioii lor the 

I'h#* ittorotnre in this 

« resu lt closer to cxpcctatlr s* 

, *irhiiiiv "on7in'-;, atwl rw'O ins o potential 

nren is vftst HIM* ro.iWly *» 

a. r,%*- nt»w oMtlnlzntlon methorts. 
sO!iroo for new o jki' '-f- 



lit 


nEFEn*=’NCES 


"“"‘"‘"Jidjiamlsch.r Auftrleb bei Ubersohrtlgesohwlndl- 

gkeit, Volta CongresSf 328-360 (1035). 

"'■'’’"swe^Siiple Conical oStrStw'Anc'' 

Lift Dependent Drag on a Slender Delta wing 

CP 428 (1^57). 
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